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SECTION I 
INTRODUCTION 



The FORTRAN Math Library consists of an extensive assortment of subroutines to aid 
the programmer in performing mathematical and trigonometric operations and functions, 
conversions between data types, bit string operations, logical relations, and other functions. 
The math routines included are for single-(real) and double -precision, complex, integer, 
and logical calculations. 

This library may be loaded in either normal or extended mode and will run in the same 
mode. 

SUBROUTINE DESCRIPTIONS 

The descriptions, in Sections IV and V, of the FORTRAN external and intrinsic functions 
and the compiler support subroutines give the name of the subroutine, its purpose, the 
DAP-16 Mod 2 calling sequence, the FORTRAN calling sequence (where appropriate), the 
method used to compute the result, the data types of the arguments and the result (where 
applicable), error messages generated by the subroutine, if any, and other routines used by 
the subroutines, if any. 

APPENDICES 

There are five appendices to this manual. Appendix A lists the contents of the library 
tapes. There are four magnetic tapes and eight paper tapes available. The first two mag- 
netic tapes are source tapes and contain the sources for the Statistical Library, FORTRAN 
Library, and Fixed Point Math Library. The third magnetic tape contains the objects for 
the software version of the three libraries; the fourth tape contains the objects for the 
hardware version of the three libraries. 



The list of "Other Routines Used" is given in the order in which they are called. If a 
routine is called more than once, it is listed only once, the first time it is called. 
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The eight object paper tapes contain the FORTRAN Library and are labeled: 

LTCF1 Tape 1 of h 

LTCF2 Tape 2 of 6 

LTCF3S Tape 3 of 6 Software Version 

LTCF3H Tape 3 of 6 Hardware Version 

LTCF4 Tape 4 of 6 

LTCF5S Tape 5 of 6 Software Version 

LTCF5H Tape 5 of 6 Hardware Version 

LTCF6 Tape 6 of 6 

The tapes are order dependent, as many of the subroutines call other subroutines which 

appear later in the library. The digits 1 through 6 on the label indicate the loading order. 

Appendix B lists the math routines by argument type. 

Appendix C lists the library subroutines by function. 

Appendix D is an alphabetical list of all the subroutines with their entry points, 
approximate storage required, subroutines referenced and the number of times referenced, 
the library tape on which they are located, and the page in this manual on which they are 
described. 

Appendix E lists the error messages produced by the subroutines and the interpretation 
of these messages. 

SYMBOLS 

The following symbols and letters are used in many of the subroutine descriptions: 

* multiplication 

/ division 

**n raised to the exponential power of n 

C complex 

D double -precision 

I integer 

L logical 

R real 

NAMING CONVENTIONS 

The intrinsic and external functions are named according to the American National 
Standards Institute (ANSI) or the Instrument Society of America (ISA) naming rules. 

The compiler support subroutines are named, for the most part, according to the 
following naming convention: The first letter of the name denotes the operation to be per- 
formed (see the list below). It is followed by a dollar sign having no significance and then 
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by two numbers. The first number (see the list below) represents the operand initially in the 
accumulator (except in load operations) and the second number represents the second operand 
or the type of result. If there is a High-Speed Arithmetic Option version of these subroutines, 
an X is appended to the name. 

Operation Argument Type 

A — Add 1 _ Integer 

C — Convert 2 — Real 

D — Divide 3 — Logical 

E — Exponential 5 — Complex 

H — Hold (store) 6 — Double -precision 

L — Load 8 — Double -precision exponent 

M — Multiply 

N — Negate 

S — Subtract 

Z — Zero (clear) 

Example s 

A$22 — add two real numbers 

D$52 — divide a complex number by a real number 

E$61 — calculate the value of a double -precision number to an 

integer power 

M$22X - multiply two real numbers, using the High-Speed 
Arithmetic Option 

LOADING INFORMATION 

There are two sets of library subroutines, one for installation with the High-Speed 
Arithmetic Option and one for those systems without this hardware option. Each set is 
contained on six rolls of paper tape. Customers who purchase the library in source form 
(on magnetic tape) receive both sets of library subroutines. 

The organization of the library is modular, thus making it possible to load only those 
routines which will be used. This concept of modularity extends to the paper tape. If 
complex or double -precision variables are not used, the first two paper tapes are not required. 

Each paper tape has been assembled via the DAP- 16 Mod 2 assembly language and 
should be loaded by the Series 16 Loader, LDR-APM. Refer to the Series 16 Equipment 
Operators' Manual, Order Number BX48, for information concerning loading object paper 
tapes. 
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SECTION II 
USE OF FORTRAN MATH LIBRARY 



DATA TYPES AND REPRESENTATIONS 

The representation of a negative number in any of the following formats (excluding logical) 
is the TWOs complement of the equivalent positive number. The complement is taken for the 
entire representation, including all subfields. The TWOs complement is taken by reversing all 
bits in the representation (ONEs complement) and adding one to the low-order position, propa- 
gating carries as required. 

Integer 

This is a 16-bit (single -precision only) word with an implied decimal point after bit 16; 
bit 1 is a sign bit (see Figure 2-1). An integer value may range from -32,768 to +32,767. 

Example: +5 = 000 000 000 000 101 = '000005 1 
-5 = 1 111 111 HI in oil = '177773 



12 3 4 5 6 7 



9 10 11 12 13 14 15 16 



V 



1 »■■-■■ I ■ » - 1 ' ■• I 



J L 



SIGN BIT 



rm 



Figure 2-1. Format of Integer 



Real 

This is a 32 -bit word in the format shown in Figure 2. 2. Bit 1 is the sign bit (0 for 
positive, 1 for negative). Bits 2-9 contain a binary number (N) with a maximum decimal value 
of 255 (377 octal) representing the 8-bit characteristic. This number is "biased" by 128 (200 
octal). The remaining 23 bits represent a binary fraction (F) with a value less than 1. The 
value represented is F*2**(N-128). A number is considered "normalized" when the fraction 
F is at least 1/2 (i. e. , the leading bit is set for a positive number). Within this representation 
the largest representable number in normalized form is just under 2**127, or approximately 
10**(38. 5). The smallest number is 2**(-129), or approximately 10**(-38. 5). The 23 magni- 
tude bits give a precision of one part in 2**23, or approximately 6. 9 digits of accuracy. Zero 



1 



The apostrophe before a number indicates octal code. 
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is shown by all zeros in these 23 bits. (Throughout this manual the word "real" is used to 
reference real singie-precision numbers.) 

Example: +5. = 100 000 111 010 000, 000 000 000 000 000 = '040720, 
-5. = 1 011 111 000 110 000, 000 000 000 000 000 = '137060, 

Double -Precis ion 

This three-word format is identical to the real number format with the exception of an 
additional 16 magnitude bits (see Figure 2-2). The 39 magnitude bits give a precision of one 
part in 2**39, or approximately 11.7 digits of accuracy. This data type should not be confused 
with hardware double-precision. 

Complex 

This is represented by two real number pairs, each having the format of a real number 
(see Figure 2-2). A real number takes two words of storage; the complex format requires 
four words. The first two words represent the real portion of the complex number, and the 
last two words represent the imaginary portion. 
Logical 

A logical value is shown as a word of all zeros for false and a value of one for true. In 
logical operations, any nonzero value is interpreted as true. The complement of a logical 
value changes it from to 1 or 1 to 0. 



FIRST WORD: 
(REAL AND 
DOUBLE- 
PRECISION) 



SECOND WORD: 
(REAL AND 
DOUBLE - 
PRECISION) 



12 3 4 



5 6 7 8 9 10 11 12 13 14 15 16 



_L__1 I I L 



II I L 



-^i^— 



CHARACTERISTIC 
(EXCESS - 128 NOTATION) 
-SIGN BIT 



MOST SIGNIFICANT 

7 BITS OF THE FRACTION 

L POSITION OF THE BINARY 
POINT 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
I ■ I 1 I I 1 I I 1 I 1 1 I I J 



v^. 



""V 



S 



N 



EXT MOST SIGNIFICANT 16 BITS OF THE FRACTION 



12 3 4 5 6 7 8 9 10 11 12 13 14 15 16 



THIRD WORD: 
(USED FOR 
DOUBLE- 
PRECISION ONLY! 



1 



1 I I L_l 1 I I I 1 1 L 



^s 



V 

LEAST SIGNIFICANT 16 BITS OF THE FRACTION 



Figure 2-2. Format of Real and Double -Precis ion Numbers 
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NORMALIZATION 

A real, double-precision, or complex number is defined as normalized when the fractional 
part has a value between 1/2 and 1. For instance, 3/8 x 2 3 and 3/4 x 2 2 both have the same 
value, but the latter is the normalized form. 

REGISTER USE 

All registers are presumed to be available to the subroutine library, and the user is 
cautioned not to expect any of them to be preserved, whether or not the arguments or results 
are stored in them. That is, any registers not specifically described as containing a particular 
result upon exit from the subroutine must be considered as having become undefined by the 
execution of the subroutine. 

ACCUMULATORS 
Integer Accumulator 

The A-register is used in all integer operations. 

Real Accumulator 

The A- and B-registers are used in all real operations. 

Complex (pseudo) Accumulator 

This four-word area in memory (AC1-AC4) is provided by the library to be used in all 
complex operations. The real portion of the complex number is stored in locations AC1 and 
AC2; the imaginary portion is stored in locations AC3 and AC4. 

Double-Precision (pseudo) Accumulator 

This three-word area in memory (AC1-AC3) is provided by the library to be used in all 
double-precision operations. 

RESULTS 

Results are stored according to their data types. Complex numbers are found in the 
complex accumulator upon exit from any of the compiler support subroutines; double -precision 
numbers are found in the double -precision accumulator; real numbers are found in the A- and 
B-registers; and integer and logical values are found in the A-register. 
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SECTION III 
DAP- 16 MOD 2 PROGRAMMING INFORMATION 



LIBRARY CALLS 

The DAP- 16 Mod 2 calling sequences for entry into the subroutines are shown in the 
descriptions in Sections IV, V and VI. When the FORTRAN compiler encounters either a 
function reference or a call to a subroutine, the following steps are performed: 

1. A call to the function or subroutine is generated. 

2. The address of each argument is determined and saved, in the order in 
which it is retrieved. In the case of expressions, this address is the 
location containing the current value of the expression. 

3. If there are two or more arguments, the final address is followed bv a 
word of zeros to serve as an argument list terminator. 

The code generated by a subprogram definition written in FORTRAN includes a call to 
the special subroutine F$AT (Argument Transfer; refer to Section VI). This call immediately 
follows the entry point and in turn is followed by a word containing a count of the number of 
arguments as defined in the definition statement, followed by that number of words. The 
F$AT subroutine fills in those words with the argument addresses (from the call to the 
subprogram) and sets the return to the word following the argument terminator word (zeros). 
All levels of indirect addressing are removed in passing these addresses. In the case of a 
single argument, the terminator word is eliminated, the argument to F$AT shows a single 
argument, and the search for the terminator is not performed. 

Null arguments may be included in a calling sequence by use of DAC*0 as the address in 
the call. Subroutines serviced by F$AT find the address DAC*0 placed in the list of addresses 
and therefore know that the parameter was null. It is equally effective to use a DAC*PTR, 
where PTR is a DAC*0. This permits a dummy argument to be null, i.e. , an argument 
passed through an intermediate subroutine call. 

The DAP- 16 Mod 2 programmer can generate his own code, performing the same functions 
as the F$AT subroutine. 

Some of the FORTRAN Math Library subroutines have additional arguments in the A- and 
B-registers, or the C-register, or the pseudo -accumulators AC1-AC4. When this is the case, 
the description references an "implicit" argument, i.e., one whose address is not explicitly 
part of the calling sequence. 

3-1 AM74 



Compiler support subroutines are those which are not normally explicitly called by the 

VDRTR AN nroarammer. For example, the statement 

Z = X + Y 

produces the following DAP- 16 Mod 2 code: 

CALL L$22 \ loadg the value of x in the A . and B-registers 

DAC X J 

A$22 \ adds the value of y to the A- and B-registers 

Y J 



CALL 

DAC 

CALL H$22 

DAC Z ' location Z 



\ stores the result in the A- and B-registers in 

J lnrafinn 7. 



Subroutines L$22, A$22, and H$22 are compiler support subroutines. They may be called 

explicitly by the FORTRAN programmer , if desired, as follows: 

CALL L$22(X) 

CALL A$22(Y) 

CALL H$22(Z) 

To perform the same function as the statement Z - X + Y and to generate the same code. 

Any of the compiler support subroutines may be called by the FORTRAN programmer in 

the following manner: 

CALL ROUTINE NAME (ARG1) 

CALL ROUTINE NAME (ARG1, ARG2) 

CALL ROUTINE NAME (ARG1, ARG2, . . . , ARGn) 

EXAMPLES OF CALLS TO LIBRARY 

CALL M$55 

DAC ARG1 

Return 
This call enters the complex multiplication subroutine, multiplying the contents of the 
complex pseudo -accumulator by the complex value in locations ARG1-ARG1+ 3 in the standard 
format for complex numbers. The result is stored in the complex accumulator (AC1-AC4), 
and any of the other registers should be presumed to have become undefined. 

CALL AMINO 

DAC I 

DAC J 

DAC K 

OCT 

Return 
This subroutine compares the three integer arguments I, J, and K (no implicit arguments) 
and returns with the value of the smallest of these, converted to data type real, in the A- and 
B-registers. Other registers are now nresumed to be undefined. 
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SECTION IV 
INTRINSIC AND EXTERNAL FUNCTIONS AND SUBROUTINES 



This section describes the mathematical and trigonometric functions and special 
FORTRAN subroutines, arranged in alphabetical order by subroutine name. 
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Purpose 



To generate the absolute value of a real number. 



DAP Calling 
Sequence 



CALL ABS 

DAC ARG1 (a real number) 
(Return) 



FORTRAN Reference ABS(R) 



Method 



This subroutine checks the real agrument, ARG1, for its algebraic 
sign. If the sign is negative, the TWOs complement of ARG1 is 
calculated. If the sign is positive, the number remains unchanged. 



Data Type of 
Arguments and 
Results 



This absolute value function of a real number gives a real 
result. 



Other Routines 
Used 



L$22, N$22 
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AIMAG 



Purpose 



To obtain the imaginary part of a complex argument and convert it 
to real format. 



DAP Calling 
Sequence 



CALL AIMAG 

DAC ARG1 (a complex number) 
(Return) 



FORTRAN Reference AIMAG (C) 



Method 



The complex argument, ARG1, is placed in the complex accumu- 
lator. The imaginary part of the complex number (AC3 and AC4) 
is then loaded into the A- and B-registers. 



Data Type of 
Arguments and 
Results 



The imaginary part of the complex argument, ARG1, is converted 
to a real number and placed in the A- and B-registers. 



Other Routines 
Used 



L$55, L$22, AC3 
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AINT 



Purpose To truncate the fractional bits of a real number. 

DAP Calling CALL AINT 

Sequence DAC ARG1 (a real number) 

(Return) 

FORTRAN Reference AINT(R) 

Method A constant (2**22) is successively added and subtracted from ARG1. 

The available precision of real numbers is such that the fractional 
part of this result is lost. If ARG1 is negative, its TWOs comple- 
ment is taken before the addition and subtraction take place and it 
is re complemented before the subroutine exits. The resultant value 
is effectively the largest integer < | ARGl| with the sign of ARG1. 



Data Type of The real argument remains a real number. 

Arguments and 

Results 

Other Routines L$22, N$22, A$22, S$22 

Used 
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ALOG 



Purpose 



To calculate the natural (base e) or common (base 10) logarithm 
of a real number. 



DAP Calling 
Sequence 



CALL ALOG 
DAC ARG 1 
(Return) 



(or ALOG10) 
(a real number) 



FORTRAN Reference ALOG(R) or ALOGIO(R) 



Method 



The log 2 of the argument, ARG1, is computed. This value is 
then converted to the desired base by multiplication by an appro- 
priate constant. 

log 2 ARG1 =F *(Cl+T(C3+T(C5+T(C7+T(C9)))))+B-.5 

where T = F *F and CI = .28853901E1 

C3 = . 96179665E0 

C5 = .57708664E0 

C7 = .41153510EO 

C9 = . 34280712E0 



" 2 

F+V^_ 
2 

F is the fractional part of the normalized argument and B is the 
binary exponent of the original argument which has been converted 
to a real number. 



Data Type of 
Arguments and 
Results 



The argument and the results are both real numbers. 



Error Messages 



The message "LG" is reported if a negative or zero-valued argu- 
ment is used, and the result is undefined. 



Other Routines 
Used 



ARG$, C$12, A$22, M$22, S$22, F$ER, H$22, L$22, D$22 
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ALOGX 



Purpose 



To calculate the natural (base e) or common (base 10) logarithm 
of a real number. 



DAP Calling 
Sequence 



CALL ALOGX (or ALOG or ALOG1 0) 
DAC ARG1 (a real number) 

(Return) 



FORTRAN Reference ALOG(R) orALOGlO(R) 



Method 



log Z = (log Z)*(log 2), where Z=ARG1. Thus for the natural 
logarithm, 

In Z = (log Z)*(log 2); for the common logarithm, 
log Z = (log Z)*(log 2). The calculation simplifies in both 
cases to a computation of log Z. Remembering that the floating- 
point number Z can be expressed as Z = F*2**B, where F is the 
fractional part and B the binary exponent of the normalized argu- 
ment Z , 



log 2 Z 



(ln(F)/ln (2))+B. 



Now, let F = F*K/K, where K may be the product 
t 

it K. 

i 

i=l 

such that F*K = 1+G,- where G is positive 

log_ X = ln(F*K/K) 

2 ln(2) + B 



ln(F*K) - ln(K) 
ln(2) 



+ B then ln(K) is 



2 ln(K ) 
i = l 





ln(2) ln(2) 


+ B 




= ln(l+G) - ln(K) 
ln(2) ln(2) 


+ B 




= G - 1/2G 2 + 1/3G 3 - 




ln(2) 




Since 


ln(2) = .69314718, 





ln(K) 
ln(2) 



+ B 



log X = 1.442695141 G -. 7213475704 G +. 4808984995 G" 



ln(K) 
ln(2) 



+ B 
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ALOGX cont. 



Data Type of 
Arguments and 
Results 



This function with a real argument results in a real number. 



Error Messages 



The message "LG" is reported if a negative or zero-valued argu- 
ment is used, and an undefined result is returned. 



Other Routines ARG$, C$12, A$22X, M$22X, S$22X, F$ER 



Used 
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ALOGIO 



Purpose To calculate the common (base 10) logarithm. 



See ALOG or ALOGX. 
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AMAXO 



Purpose To find the maximum real value in a list of intege 



rs. 



See MAXO. 
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AMAX1 



Purpose To find the maximum real value in a list of real arguments. 



See MAX1. 
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AMINO 



Purpose To find the minimum real value in a list of integers. 



See MINO. 
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AMIN1 



Purpose To find the minimum real value in a list of real arguments. 



See MINI. 
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AMOD 



Purpose 



To compute the remainder resulting from the division of two real 
numbers. 



DAP Calling 
Sequence 



CALL AMOD 
DAC ARG1 
DAC ARG2 
OCT 
(Return) 



(real dividend) 

(real divisor) 

(end of arguments flag) 



FORTRAN Referenc e 



AMOD(R, R) 



Method 



This subroutine divides ARG1 by ARG2 by calling D$22. The functior. 
AMOD (ARG1, ARG2) is defined as: 

Al - (A1/A2) *A2, where A1=ARG1 and A2=ARG2 

(A1/A2) is the integer whose magnitude does not exceed the magni- 
tude of A1/A2 and whose sign is the same as that of A1/A2. 



Data Type of 
Arguments and 
Results 



This function with two real arguments results in a real number for 
a remainder. 



Other Routines 
Used 



L$22, D$22, AINT, M$22, N$22, A$22 
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ATAN 



Purpose 



To calculate the principal value of the arctangent (i. e. , 1st or 4th 
quadrant angle) of a real number or to compute and adjust for 
quadrant the arctangent of a real number expressed as a ratio 
(X/Y). 



DAP Calling 
Sequence 



CALL ATAN or CALL ATAN2 

DAC ARG1 (a real number) DAC ARG1 (both arguments 

(Return) DAC ARG2 are real numbers 

OCT (end of arguments flag) 
( Return) 



FORTRAN Reference ATAN(R) or ATAN2(R, R) 



Method 



For ATAN, let N = ABS(ARGl).. The arctangent of N is evaluated by- 
dividing the total range 0<N<10**75 into three intervals: 

If N<10**(-B), ATAN(N) =N 

If N>10**10, ATAN(N) = pi/2 

If 10** (-8) < N < 10**10, 

ATAN(N) = base angle + P(Z) 

= base angle + C1*Z+C3*Z**3+C5*Z**5. . . 

If N < 1 /2, Z - N and base angle * 

If N<2, Z = (N-1)/(N+1) and base angle = pi/4 

If N<2, Z = (-1/N) and base angle = pi/2 

For ATAN2, the arctangent of the quotient of ARG1/ARG2 (ARG1 = 
side opposite, ARG2 = side adjacent, or sin/cos) is computed as 
in ATAN and adjusted for quadrant by examination of the signs of 
the numerator and denominator. 



Results 



Quadrant ARG1 ARG2 Quotient 



Results (radians) 



1 


+ 


+ 


to 00 


to pi/2 


2 


+ 


- 


- to 


pi/2 to pi 


3 


- 


- 


to 00 


-pi to -pi/2 


4 


. 


+ 


- to 


-pi/2 to 
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ATAN cont. 



Data Type 


of 


Arguments 


and 


Results 




Other Routines 



Used 



This arctangent function of a real number results in a real numbe 



ARG$, D$22, N$22, M$22, A$22, S$22 
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ATAN2 



Purpose To calculate the arctangent as the quotient of two real numbers. 



See ATAN. 



4-16 A M74 



Purpose 



To generate the absolute value of a complex number. 



DAP Calling 
Sequence 



CALL CABS 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CABS(C) 



Method 



The argument is squared and its square root is taken to arrive at 
its absolute value; e. g. , if ARG1 = X+IY, 
CABS(ARGl) = SQRT (X**2+Y**2). 



Data Type of 
Arguments and 
Results 



This absolute value function of a complex number gives a real 
result. 



Other Routines F$AT, SUB$, L$22, M$22, H$22, A$22, SQRT 



Used 
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Purpose 



To calculate the cosine of a complex number with the real part 
in radian measure. 



DAP Calling 
Sequence 



CALL CCOS 
DAC ARG 1 

(Return) 



(a complex number) 



FORTRAN Reference 



CCOS(C) 



Method 



The cosine function is transtormed into the sine function by use of 
the trigonometric identity COS (Z) = SIN (Z + pi/2), where Z=Y+IY. 
SIN (Z+pi/2) is then evaluated. 



Data Type of 
Arguments and 
Results 



This cosine function of a complex number results in a complex number. 



Other Routines 
Used 



F$AT,L$55, A$55, H$55, CSIN 



4-18 



AM74 



Purpose 



To calculate the exponential of a complex number with the imaginary 
part in radian measure. 



DAP Calling 
Sequence 



CALL CEXP 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CEXP(C) 



Method 



The following algorithm is used to calculate the value of e**ARGl, 
where ARG1 is a complex number: 

If ARG1 = X+IY, 

e**(X+IY) = (e**X) * (e**IY) = (e**X) *COS(Y) + I * (e**X) *SIN(Y) 



Data Type of 
Arguments and 
Results 



This function raises e to a complex power and gives a complex result. 



Other Routines F$AT, SUB$, EXP, H$22, COS, M$22, SIN, L$55 

Used 
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CLOG 



Purpose 



To calculate a particular value of the natural logarithm (base e) 
of a complex number. 



DAP Calling 
Sequence 



CALL CLOG 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CLOG(C) 



Method 



The following algorithm is used to calculate ln(ARGl), where 
ARG1 = X + IY: 

In 'X 4 - IY) - P + T( * ) 

where R = In (X**2+ Y**2)**. 5 = 1/2 In (X**2+Y**2) 

4> = (TAN**-1)(Y/X) = 4> + 2Kpi 

where K = 0, ± 1, ± 2, ... 

A particular value for <p is chosen such that -pi<^>< pi by enter- 
ing the arctangent routine ATAN2 . 



Data Type of 
Arguments and 
Results 



This logarithm function of a complex number gives a complex result. 



Other Routines 
Used 



F$AT, L$22, M$22, H$22, A$22, ALOG, ATAN2, L$55, SUB< 
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CMPLX 



Purpose 



To combine two real numbers into one complex quantity. 



DAP Calling 
Sequence 



CALL CMPLX 

DAC ARG1 (a real number) 

DAC ARG2 (a real number) 

OCT (end of arguments flag) 

(Return) 



FORTRAN Reference CMPLX (R, R) 



Method 



The first real argument (ARG1) is stored in the real portion of the 
complex accumulator (AC 1 and AC2). The second real argument 
(ARG2) is stored in the complex portion of the complex accumulator 
(AC3 and AC4). 



Data Type of 
Arguments and 
Results 



The two real arguments are combined into one complex number and 
stored in the complex accumulator. 



Other Routines 
Used 



F$AT, SUB$, L$22, H$22, L$55 
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CONJG 



Purpose 



To obtain the conjugate of a complex number. 



DAP Calling 
Sequence 



CALL CONJG 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CONJG(C) 



Method 



This subroutine reverses the sign of the imaginary part of the 
complex argument (ARG1). 



Data Type of 
Arguments and 
Results 



The complex argument in this function remains a complex number. 



Other Routines F$AT, SUB$, L$22, H$22, N$22, L$55 

Used 
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Purpose To calculate the cosine of a real number expressed in radians. 



See SIN. 
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CSIN 



Purpose 



To calculate the sine of a complex number with the real part in 
radian measure. 



DAP Calling 
Sequence 



CALL CSIN 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CS1N(C) 



Method 



The sine function of the complex number ARG1 (X+IY) is computed 
as follows: 

SIN (X+IY) = SIN(X) * COSH (Y) + I * (COS(X) * SINH(Y) ) 

where SINH(Y) = 1/2 * (E**Y-E**-Y) 

COSH(Y)= 1/2 * (E**Y+E**-Y) 



Data Type of 
Arguments and 
Results 



The argument and the result of this function are complex numbers. 



Other Routines 
Used 



F$AT, SUB$, EXP, H$22, L$22, D$22, A$22, SIN, M$22, S$22, COS, 
L$55 
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Purpose 



To calculate the square root of a complex number. 



DAP Calling 
Sequence 



CALL CSORT 

DAC ARG1 (a complex number) 

(Return) 



FORTRAN Reference CSQRT(C) 



Method 



If the complex argument is positive, (A+B)**. 5 - C+DI is determined 
as follows: 

C = (((A**2+B**2)**.5+A)/2)**. 5 

D = B/(2*C) 

If the argument is negative, ABS(D) = ( ( (A**2+B**2) -A)/2)**. 5. 

The sign of the real part of the result will be positive and the sign 
of the imaginary part of the result will be the same as the sign of 
the imaginary part of the argument. That is, the results will lie 
in quadrants I or IV of the complex plane. 



Data Type of 
Arguments and 
Results 



This square root function of a complex number results in a complex 
number. 



Other Routines 
Used 



F$AT, SUB$, CABS, H$22, ABS, A$22, M$22, SQRT, L$22, D$22, 
L$55 
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Purpose To generate the absolute value of a double -precision number. 



DAP Calling CALL DABS 

Sequence DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference DABS(D) 



Method This subroutine checks the double -precision argument, ARG1, for 

its algebraic sign. If the sign is negative, the TWOs complement o 
ARG1 is calculated. If the sign is positive, the number remains 
unchanged. 



Data Type of This function with a double -precision argument results in a 

Arguments and double -precision number. 

Results 



Other Routines F$AT, L$66, N$66 

Used 
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DATAN 



Purpose 



To calculate the arctangent of a double -precision number. 



DAP Calling 
Sequence 



CALL DATAN 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference DATAN (D) 



Method 



The principal value is computed. See "Method" for ATAN. 



Data Type of 

Arguments and 
Results 



This function with a double -precision argument results in a double - 
precision number. 



Other Routines F$AT f DABS, H$66, C$81, L$66, A$66, D$66, M$66, N$66 



Used 
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DATAN2 



Purpose 



To calculate the arctangent of the quotient of two double -precision 
numbers. 



DAP Calling 
Sequence 



CALL DATAN2 

DAC ARG1 (a double -precision number (X) ) 

DAC ARG2 (a double -precision number (Y) ) 

OCT (end of arguments flag) 

(Return) 



FORTRAN Reference DATAN2(D,D) 



Method 



The arctangent of the quotient (X/Y) is adjusted for the quadrant 
by examining the signs of the numerator and denominator. See 
"Method" for ATAN. 



Data Type of 
Arguments and 
Results 



This arctangent function of a double -precision quantity gives a 
double -precision result. 



Error Messages 



The error message "DT" is reported if the second argument is 
zero. The result in the double -precision accumulator is undefined. 



Other Routines 
Used 



F$AT, L$66, H$66, F$ER, D$66, DATAN, S$66, A$66 
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Purpose 



To convert a real number to double -precision format. 



DAP Calling 
Sequence 



CALL DBLE 

DAC ARG1 (a real number) 

(Return) 



FORTRAN Reference DBLE(R) 



Method 



This subroutine stores the real argument, ARG1, in AC1 and AC2. 
A word of zeros is appended to the real number as the least signifi- 
cant word of the double -precision fraction and stored in AC 3. 



Data Type of 
Arguments and 
Results 



The real argument is converted to a double -precision number. 



Other Routines F$AT, L$22, C$26 

Used 
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Purpose To calculate the cosine of a double -precision number expressed in 

radians. 



DAP Calling CALL DCOS 

Sequence DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference DCOS(D) 



Method The cosine function is transformed into the sine function using the 

trigonometric identity COS (X) = SIN (pi/2+X). SIN (pi/2+X) is then 
evaluated, with X = ARG1. 



Data Type of This function with a double -precision argument gives a double - 

Arguments and precision result. 

Results 



Other Routines F$AT, L$66, A$66, H$66, DS1N 

Used 



4-30 A.M7<= 



Purpose 



To calculate e**x, where x is a double-precision number. 



DAP Calling 
Sequence 



CALL DEXP 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference DEXP(D) 



Method 



In calculating e**ARGl, the following method is used- e**ARGl 
2**(ARGl*log 2 (e)) = 2**(I+F), where land F are the integer and 
fractional portions, respectively, of the product ARGl*log 2 (e). 



Data Type of 
Arguments and 
Results 



This function raises e to the power of a double-precision argument 
and gives a double-precision result. 



Other Routines 
Used 



F$AT, L$66, M$66, H$66, C$61, C$16, N$66, A$66, S$66, D$66 
A$81 ' 
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DIM 



Purpose 



To compute the positive difference between two real arguments. 



DAP Calling 
Sequence 



CALL DIM 
DAC ARG1 
DAC ARG2 
OCT 

(Return) 



(a real number) 
(a real number) 
(end of arguments flag) 



FORTRAN Reference DIM(R, R) 



Method 



ARG1 - ARG2 is computed. If the result is positive, this value is 
the result given. If ARG1 - ARG2 is a negative quantity, the 
result of this function is zero. 



Data Type of 
Arguments and 
Results 



This routine to calculate the difference between two real numbers 
results in a real number. 



Other Routines 
Used 



L$22, S$22 



4-32 



AM74 



DINT 



Purpose 



To truncate the fractional bits of a double-precision number. 



DAP Calling 
Sequence 



CALL DINT 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference 



DINT(D) 



Method 



A constant (2**38) is successively added and subtracted from the 
argument, ARG1. The available precision of double -precision 
numbers (39 bits) is such that the fractional part of this result is 
lost. If ARG1 is negative, its TWOs complement is taken before 
the addition and subtraction take place and it is recomplemented 
before the subroutine exits. The resultant value is effectively 
the largest integer <|ARG1 | with the sign of ARG1. 



Data Type of 
Arguments and 
Results 



The double -precision argument after truncation remains a double- 
precision number. 



Other Routines 
Used 



L$66, N$66, A$66, S$66, AC1 
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DLOG 



Purpose 



To calculate the natural (base e) logarithm of a double-precision 

number, 



DAP Calling 
Sequence 



CALL DLOG 
DAC ARG1 

(Return) 



(a double -precision number) 



FORTRAN Reference 



DLOG(D) 



Method 



This routine is also used by DLOG2 and DLOG10. Log A (X), where 
X = ARG1, is calculated as log 2 (X)/log 2 (A). To calculate log 2 (X), 
X is considered as the number F 1 *(2**B) where 1/2 <F < 1. Log 2 
(X) = log 2 (F 1 ) + the binary exponent of F , and log 2 (F ) = 1/2 + 
C1*Z + C3(Z**3) + ... where 

,1 

CI = 2.885390081845024D0 

C3 = .9617966484737566D0 

C5 = . 577086624639535D0 

C7 = .4115350984570017D0 

C9 = .3428071228932386D0 



-Sz) 

(Fl + Jl) 



Data Type of 
Arguments and 
Results 



This natural logarithm function of a double -precision argument 
results in a double-precision number. 



Error Messages 



The message "DL" is reported if a negative or zero -valued argu- 
ment is found. The result in the double -precision accumulator is 
undefined. 



Other Routines 
Used 



F$AT, DLOG2, M$66 
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Purpose 



To calculate the common (base 2) logarithm of a double -precision 

number. 



DAP Calling 
Sequence 



CALL DLOG2 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference DLOG2(D) 



Method 



This routine is used by DLOG and DLOG10 to calculate log ? (X), 
where X is equal to F**(2**B) and 1/2 sF< 1. See "Method" for 
DLOG. 



Data Type of 
Arguments and 
Results 



This common logarithm function with a double -precision argument 
results in a double-precision number. 



Error Messages 



The message "DL" is reported if a negative or zero-valued argu- 
ment is found. The result is undefined. 



Other Routines 
Used 



F$AT, L$66, F$ER, C$81, C$16, H$66, Z$80, A$66, S$66, D$66, 
M$66 
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DLOGIO 



Pu rpo s e 



To calculate the common (base 10) logarithm of a double-precision 
number. 



DAP Calling 
Sequence 



CALL DLOGIO 

DAC ARG1 (a double-precision number) 

(Return) 



FORTRAN Reference DLOGIO(D) 



Method 



See "Method" for DLOG. 



Data Type of 
Arguments and 



This logarithm function with a double -precision argument results 
in a double -precision number. 



Error Messages 



The message "DL" is reported if a negative or zero-valued argu- 
ment is found. The result is undefined. 



Other Routines 
Used 



F$AT, DLOG2, M$66 
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DMAX1 



Purpos 



To find the largest value in a list of double-precision arguments. 



DAP Calling 
Sequence 



CALL DMAX1 

DAC ARG1 (first double -precision argument) 

DAC ARG2 (a double -precision number) 



DAC ARGn (last double -precision argument) 

OCT (end of arguments flag) 

(Return) 



FORTRAN Reference DMAX1 (D,D D) 



Method 



Compare the arguments and retain the largest value. 



Data Type of 
Arguments and 
Results 



The largest double-precision argument is stored in the double- 
precision accumulator. 



Other Routines 
Used 



L$66, H$66, S$66 
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DMIN1 



Purpose 



To find the smallest value in a list of double-precision arguments. 



DAP Calling 
Sequence 



CALL DMIN1 
DAC ARG1 

DAC ARG2 



(a double -precision argument) 
(a double -precision argument) 



DAC ARGn 

OCT 

(Return) 



(last double-precision argument) 
(end of arguments flag) 



FORTRAN Reference DMIN1 (D,D, . . . ,D) 



Method 



Compare the arguments and retain the smallest value. 



Data Type of 
Arguments and 
Results 



Both of the arguments are double -precision and the result of this 
function is a double-precision number. 



Other Routines 
Used 



L$66, H$66, S$66 



4-38 



AM 74 



DMOD 



Pu rpo s e 



To compute the remainder resulting from the division of two double- 
precision numbers. 



DAP Calling 
Sequence 



CALL DMOD 
DAC ARG1 
DAC ARG2 
OCT 
(Return) 



(a double precision number) 
(a double-precision number) 
(end of arguments flag) 



FORTRAN Reference 



DMOD(D, D) 



Method 



This subroutine divides ARG1 by ARG2 by calling D$66. The 
function DMOD (A1.A2) is defined as Al -(A 1/A2)*A2, where (A1.A2) 
is the integer whose magnitude does not exceed the magnitude of 
A1/A2 and whose sign is the same as that of A1/A2. 



Data Type of 
Arguments and 
Results 



This function with two double-precision arguments results in a 
double-precision number for a remainder. 



Other Routines 
Used 



F$AT, L$66, D$66, H$66, DINT, M$66, S$66, N$66 



4-39 



AM 74 



DSIGN 



Purpose 



To generate a value consisting of the sign of the second double- 
precision argument and the magnitude of the first double-precision 
argument. 



DAP Calling 
Sequence 



CALL DSIGN 
DAC ARG1 

DAC ARG2 

OCT 

(Return) 



(a double-precision number) 
(a double -precision number) 
(end of arguments flag) 



FORTRAN Reference 



DSIGN(D, D) 



Method 



ARG2 is tested for its algebraic sign and, depending on the sign of 
ARG1, the procedure is as follows: 



ARG1 


ARG2 


Result 


_ 


+ 


+ |ARG1 


_ 


- 


- |ARG1 


+ 


+ 


+ |ARG1 


+ 


_ 


- IARG1 



Data Type of 
Arguments and 
Results 



Both arguments for this call are double-precision numbers and the 
result is a double -precision number. 



Other Routines 
Used 



F$AT, L$66, N$66 
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DSIN 



Purpose 



To calculate the sine of a double -precision number expressed in 
radians. 



DAP Calling 
Sequence 



CALL DSIN 
DAC ARG1 

(Return) 



(a double -precision number) 



FORTRAN Reference 



DSIN(D) 



Method 



An arbitrary angle X expressed in radian measure can be reduced 
to the range < Y < jp- through the relation X = Y + N(pi/2). 
Adjustment is made for quadrant before using a modified Taylor's 
expansion. 



Data Type of 
Arguments and 
Results 



This sine function with a double -precision argument results in a 
double -precision number. 



Other Routines 
Used 



F$AT, DABS, M$66, H$66, C$61, C$16, N$66, A$66, MOD, 
L$66, S$66 
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Purpose 



To calculate the square root of a double -precision number. 



DAP Calling 
Sequence 



CALL DSQRT 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference 



DSQRT (D) 



Method 



A first approximation to the double -precision square root of the 
double -precision argument is obtained by calling the real square 
root routine (SQRT). One more Newton -Raphs on iteration is the 
made to achieve full double -precision accuracy. 



Data Type of 
Arguments and 
Results 



This square root function of a double -precision argument result 
in a double -precision number. 



Other Routines 



F$AT, L$66, C$62, H$22, SQRT, C$26, H$66, D$66, A$66, 
A$81 
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Purpose 



To calculate e**x, where x is a real number. 



DAP Calling 
Sequence 



CALL EXP 

DAC ARG1 (a real number) 

(Return) 



FORTRAN Reference 



EXP(R) 



Method 



e**ARGl = 2**(ARG1 * log^e)) = 2**(I+F), where I is the integer 
and F is the fractional portion of the product ARG1 log (e). The 
value of F is used to define the quantities I 1 , F(l), and F(2): 



F(i; 



-1 < F < -1/2 


I - I 


1/4 


-1/2 < F < 


I - I 


3/4 


< F < 1/2 


I 


1/4 


1/2 <F < 1 


I 


3/4 



F(2) 
F + 3/4 
F + 1/4 
F - 1/4 
F - 3/4 



From the above table, e**ARGl=2**(I'+Fl+F2) = 2**(I'+F1) * (2**F2) 

where 

2**F2 = e**(F2*ln(2)) = e**F = (A(F) )/(A(F)-B(F) ) 
A(F) = CH(F*F), B(F) = C2*F 



Data Type of 
Arguments and 
Results 



id 
This exponential function with a real argument (e ) results in a real 
number. 



Error Message 



When overflow occurs, the error message "EX" is reported and the 
answer returned is the maximum value possible (1. 7E38). When 
underflow occurs, the value is returned without an error message. 



Other Routines 
Used 



ARG$, N$22, M$22, S$22, A$22, D$22, F$ER 
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Purpose 



To convert an integer argument to real format. 



DAP Calling 
Sequence 



CALL FLOAT 

DAC ARG1 (an integer value) 

(Return) 



FORTRAN Reference FLOAT(I) 



Method 



This routine extracts the integer and converts it to real format, 
leaving the result in the A- and B-registers. 



Data Type of 
Arguments and 
Results 



This routine converts an integer argument to a real number. 



Other Routines 
Used 



C$12 
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IABS 



Purpose To generate the absolute value of an integer. 



DAP Calling CALL IABS 

Sequence DAC ARG1 (an integer value) 

(Return) 



FORTRAN Reference IABS (I) 



Method 



This subroutine checks the integer argument, ARG1, for its algebraic 
sign. If the sign is negative, the TWOs complement of ARG1 is cal- 
culated. If the sign is positive, the number remains unchanged. 



Data Type of This absolute value function with an integer argument results in an 

Arguments and integer. 

Results 
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IDIMI 



Purpose 



To compute the positive difference between two integer arguments. 



DAP Calling 
Sequence 



CALL IDIM 
DAC ARG 1 
DAC ARG 2 
OCT 
(Return) 



(an integer value) 
(an integer value) 
(end of arguments flag 



FORTRAN Reference IDIM(I, I) 



Method 



Compute DIF = ARG1-ARG2. If DIF is positive, the result of this 
function is the value of DIF. If DIF is negative, the result of this 
function is zero. 

DIF = ARG1 - MIN(ARG1, ARG2) 



Data Type of 
Arguments and 
Results 



The result of this function with two integer arguments is an integer. 
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IDINT 



Purpose To truncate the fractional bits from a double-precision argument, 

thus converting it to integer format. 



See IFIX. 



4 " 47 AM74 



Purpose 



To fetch the contents of the memory location specified by ARG1. 



DAP Calling 
Sequence 



CALL 1FETCH 
DAC ARG 1 

(Return) 



FORTRAN Reference IFETCH(ARGl) 



Method 



The A-register is loaded with the contents of the location specified 
by ARG 1 . 



Other Routines ARG$ 

Used 
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IFIX 



Purpose 



To truncate the fractional bits from a real or double -precision 
argument, thus converting it to integer format. 



DAP Calling 
Sequence 



CALL IFIX (or CALL INT) 

DAC ARG1 (a real number) 

(Return) 



CALL IDINT 

DAC ARG1 (a double -precision number) 

(Return) 



FORTRAN Reference 



IFIX(R), INT(R), IDINT (D) 



Method 



This subroutine truncates the fractional bits of ARG1, shifts it to 
the right until the binary point is at the end of the register, and 
normalizes the result. It then uses the characteristic to scale 
the value to an integer. 



Data Type of 
Arguments and 
Results 



If either IFIX or INT is called, the argument is a real number and 
the result is an integer. If IDINT is called, the argument is a double- 
precision number and the result is an integer. 



Other Routines 
Used 



L$22, C$21 
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INT 



Purpose To truncate the fractional bits from a real argument, thus convert- 

ing it to integer format. 



See IFIX. 
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ISIGN 



Purpose 



To generate a value consisting of the sign of the second integer 
argument and the magnitude of the first integer argument. 



DAP Calling 
Sequence 



CALL ISIGN 
DAC ARG1 

DAC ARG2 

OCT 

(Return) 



(an integer value) 
(an integer value) 
(end of arguments flag) 



FORTRAN Reference 



ISIGN(I, I) 



Method 



ARG2 is tested for its algebraic sign and, depending on the sign of 
ARG1, the procedure is as follows: 



ARG1 



ARG2 



Result 



+ 


+ 


+ IARG1 1 


+ 


- 


- |ARG1 | 


- 


+ 


+ IARG1 | 


- 


- 


- IARG1 | 



Data Type of 
Arguments and 
Results 



Both arguments and the result are integers. 
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Purpose To store the contents of the second argument in the location speci- 

fied as the first argument. 



DAP Calling CALL ISTORE 

Sequence DAC ARG1 (target word address) 

~ DAC ARG2 (word to be stored) 

OCT (end of arguments flag) 

(Return) 



FORTRAN Reference ISTORE(ARGl , ARG2) 



Method Fetch the target word address (ARG1) and save it. 

Fetch the word to be stored (ARG2) and use it to replace the 
contents of the target location. Effectively, the contents of ARG2 
are stored in location ARG1. 



Other Routines F$AT 

Used 
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LOC 



Purpose To determine the address of the argument. 



DAP Calling CALL LOC 

Sequence DAC ARG 1 

(Return) 



FORTRAN Reference LOC (ARG 1 ) 



Method Fetch the argument address (direct or indirect) and load it into 

the A-register. 
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MAXO 



Purpose 



To find the largest value in a list of integer arguments and exit 
with this value or convert it to real format (AMAXO) and exit. 



DAP Calling 
Sequence 



CALL MAXO (or AMAXO) 

DAC ARG1 (integer value) 

DAC ARG2 (integer value) 



DAC ARGn 

OCT 

(Return) 



(last integer argument) 
(end of arguments flag) 



FORTRAN Reference 



MAX0(I,I, ... ,1) 



or AMAX0(I,I, . . . , I) 



Method 



This subroutine compares the arguments and retains the largest 
value. If AMAXO is called, the result is converted to real by 
calling FLOAT before the subroutine exits. 



Data Type of 
Arguments and 
Results 



The arguments are integers in either call (MAXO or AMAXO). The 
result is integer if MAXO is called; the result is a real number if 
AMAXO is called. 



Other Routines 
Used 



FLOAT 
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MAX1 



Purpose 



To find the largest value in a list of real arguments and exit with 
this value or convert it to an integer (MAX1) and exit. 



DAP Calling 
Sequence 



CALL MAX1 (orAMAXl) 
DAC ARG1 (a real number) 

DAC ARG2 (a real number) 



DAC ARGn 

OCT 

(Return) 



(last real argument) 
(end of arguments flag) 



FORTRAN Reference 



MAX1(R,R,...,R) or AMAX1 (R , R, . . . , R) 



Method 



This subroutine compares the arguments and retains the largest 
value. If MAX1 is called, the result is converted to integer by 
calling IFIX before the subroutine exits. 



Data Type of 
Arguments and 
Results 



The arguments are real numbers in either call (AMAX1 or MAX1). 
The result is real if AMAX1 is called; the result is an integer if 
MAX1 is called. 



Other Routines 
Used 



L$22, H$22, S$22, IFIX 
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MINO 



Purpose 



To find the smallest value in a given set of integers and exit with 
this value or convert this value to a real number and exit. 



DAP Calling 
Sequence 



CALL MINO 
DAC ARG 1 
DAC ARG 2 



(or AMINO) 

(an integer value) 

(an integer value) 



DAC ARGn 

OCT 

(Return) 



(last integer argument) 
(end of arguments flag) 



FORTRAN Reference 



MIN0(I,I, 



,1) or AMIN0(I,I, .. .,!) 



Method 



This subroutine compares the arguments and retains the smallest 
value. If AMINO is called, the result is converted to a real number 
before the subroutine exits. 



Data Types of 
Arguments and 
Results 



The arguments are integers in either call (MINO or AMINO). The 
result is integer if MINO is called; the result is a real number if 
AMINO is called. 



Other Routines 
Used 



FLOAT 
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MINI 



Purpose 



To find the smallest value in a list of real arguments and exit with 
this value (AMIN1) or convert it to an integer (MINI) and exit. 



DAP Calling 
Sequence 



CALL MINI 
DAC ARG1 

DAC ARG2 



(or AMIN1) 

(a real number) 

(a real number) 



DAC ARGn 

OCT 

(Return) 



(last real argument) 
(end of arguments flag) 



FORTRAN Reference MIN1(R,R R) or AMIN1 (R, R, . . . , 



R) 



Method 



Compare the arguments and retain the smallest value 



Data Type of 
Arguments and 
Results 



The arguments are real numbers for either call (MINI or AMIN1). 
The result is real if AMIN1 is called; the result is integer if MINI 
is called. 



Other Routines 
Used 



L$22, H$22, S$22, IFIX 
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MOD 



Purpose To compute the remainder resulting from the division of two integers. 



DAP Calling CALL MOD 

Sequence DAC ARG1 (an integer value) 

" DAC ARG2 (an integer value) 

OCT (end of arguments flag) 

(Return) 



FORTRAN Reference MOD(I, I) 



Method This subroutine divides ARG1 by ARG2 by calling D$ll. The function 

MOD(Al,A2) is defined as Al -(Al /A2)*A2, where (A1/A2) is the 
integer whose magnitude does not exceed the magnitude of A1/A2 
and whose sign is the same as that of A1/A2. 



Data Type of This function with two integer arguments results in an integer for a 

Arguments and remainder. 

Results 



Other Routines D$ll, M$ll 

Used 
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Purpose To check for an error condition. 



DAP Calling CALL OVERFL 

Sequence DAC J (an integer value) 

(Return) 



FORTRAN Reference OVERFL(J) 



Method This subroutine checks error flag AC5 for a nonzero value, which 

indicates that an entry to the error subroutine, F$ER, was made 
since the last call to OVERFL. If AC5 is nonzero, the variable J 
is set to 1 and AC5 is cleared. If AC5 is zero, J is set to 2. 



Other Routines AC 5 

Used 
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Purpose 



To load the real portion of a complex number into the A- and re- 
register. 



DAP Calling 
Sequence 



CALL 
DAC 

(Return) 



REAL 
ARG1 



(a complex number) 



FORTRAN Reference REAL(C) 



Method 



This subroutine calls ARG$ to place the complex argument, ARG1, 
into the index register. The real portion, i.e. , the first two 
words, of the complex argument is then loaded into the A- and B- 

reeisters . 



Data Type of 
Arguments and 
Results 



This function of a complex number results in a real number. 



Other Routines 
Used 



ARG$ 
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SIGN 



Purpose 



To generate a value consisting of the sign of the second real argu- 
ment and the magnitude of the first real argument. 



DAP Calling 
Sequence 



CALL SIGN 
DAC ARG1 

DAC ARG2 

OCT 

(Return) 



(a real number) 
(a real number) 
(end of arguments flag) 



FORTRAN Reference SIGN(R, R) 



Method 



ARG2 is tested for its algebraic sign and, depending on the sign of 
ARG1, the procedure is as follows: 



ARGI 

+ 
+ 



ARG2 



Result 

+ | ARGI | 

| ARGI | 

+ | ARGI I 

| ARGI | 



Data Type of 
Arguments and 
Results 



Both arguments are real numbers and the result is a real number. 



Other Routines 
Used 



L$22, N$22 
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SIN 



Purpose 



To calculate the sine or cosine of a real number expressed in radians. 



DAP Calling 
Sequence 



CALL SIN (or COS) 

DAC ARG1 (a real number) 

(Return) 



FORTRAN Reference 



SIN(R) or COS(R) 



Method 



The angle is reduced to the first quadrant by the use of the relation 
X = Y+N*(pi/2) and the identities SIN(Y) = COS(pi/2-Y) and COS(Y) 
= SIN(pi/2-Y). A modified Taylor's expansion is then used to cal- 
culate the sine of the first quadrant angle. 

The cosine function is transformed into the sine function by the use 
of the identity COS(X) - SIN(pi/2-X); SIN(pi/2-X) is then evaluated, 
wnere X =ARG1 



Data Type of 
Arguments and 
Results 



This sine function with a r^ea. 1 argument results in a real number. 



Other Routines 
Used 



ARG$, N$22, M$22, S$22, A$22 
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Purpose 



To set or reset the pseudo sense lights and switches. 



DAP Calling 
Sequence 



CALL SLITE 

DAC ARG 1 

(Return) 



CALL 


SLITET 


DAC 


ARG1 


DAC 


ARG2 


OCT 






(Return) 



(where ARG1 is the address of the variable con- 
taining the sense light number). 

(or CALL SSWTCH) 

(where ARG1 is the address of the variable con- 
taining the sense light or switch (SSWTCH) 
number to be interrogated, and ARG2 is the 
address of the location in which to store the 
"set or reset" indicator; (l=set, 2 = reset). 



FORTRAN Reference 



CALL SLITE (I), CALL SLITET(I, J), CALL SSWTCH(I, J) 



Method 



SLITE — The ARG$ routine is used to place the variable address in 
the index register. The argument (I) is tested for zero. If zero, all 
sense light positions are reset; otherwise, the sense light specified 
is shifted to its appropriate position and INCLUSIVELY ORed with 
current settings, leaving them undisturbed. 

SLITET -The ARG$ routine is used to place the sense light number 
in the A-register and the location of the variable in the index register. 
If the sense light number is 0, a 2 is inserted into the variable J, 
signifying a reset condition. Otherwise, the sense light bit is moved 
to its proper position in the A-register. A logical AND is executed 
with the sense light register. If the result of the AND is zero, the 
sense light is reset and a 2 is placed in J. If the result of the AND 
is not zero, an EXCLUSIVE OR is carried out with the sense light 
register, resetting the sense light specified and storing a 1 in J to 
signify that the sense light was set on entry. 

SSWTCH -The ARG$ routine is used to place the sense switch number 
in the A-register and the variable location in the index register. If 
the sense switch number is (no real switch), J is set to 1 . If the 
sense switch number is valid (1 to 4), J is set to 1 if the external 
switch is set and set to 2 if the external switch is not set. 



Other Routines 
Used 



ARG$, L$33 
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Purpose To set or reset the pseudo sense lights and switches. 



See SLITE. 
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Purpose 



To calculate the square root of a real number. (This subroutine has 
a high-speed version, SQRTX. ) 



DAP Calling 
Sequence 



CALL SORT 

DAC ARG 1 

(Return) 



(a real number) 



FORTRAN Refe r enc e 



SQRT(R) 



Method 



Given the argument N = F(2**e), the mantissa is adjusted so that 
e is even and 1 /4 < e < 1 . An initial approximation to the square 
root (Y) is chosen as follows: 

Y = 7/8(F) + 9/32 if e < 1 /2 

Y = 9/16(F) + 7/16 if e > 1/2 

Two Newton-Raphson iterations are then made to obtain full single- 
precision accuracy. 



Data Type of 
Arguments and 
Results 



This square root function of a real number results in a real numbt 



Error Messages 



The error message "SQ" is reported if a negative argument is found. 
An undefined result is returned in the A-and B-registers. 



Other Routines 
Used 



ARG$, DIV$, D$22, A$22, F$ER 
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Purpose 



To calculate the square root of a real number. (This routine re- 
quires the High-Speed Arithmetic Option. ) 



DAP Calling 
Sequence 



CALL SQRTX 
DAC ARG1 

(Return) 



(or SQRT) 

(a real number) 



FORTRAN Reference 



SQRT(R) 



Method 



Given the argument N=F*(2**e), the mantissa is adjusted so that 
e is even and 1 /4 < e < 1. An initial approximation to the square 
root of ARG1 is chosen as follows: 

ARG1 = 7/8(F) + 9/32 if e < 1 /2 

ARG1 = 9/16(F) + 7/16 if e > 1/2 

Two Newton-Raphson iterations are then made to obtain full single- 
precision accuracy. 



Data Type of 
Arguments and 
Results 



This square root function of a real number results in a real numbe; 



Error Messages 



The error message "SQ" is reported if a negative argument is found. 
An undefined result is returned in the A-and B-registers. 



Other Routines 
Used 



ARG$, D$22X, A$22X, F$ER 
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Purpose To set or reset the pseudo sense switches. 



See SLITE. 
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I Jmt. Ni H 



Purpose 



To calculate the hyperbolic tangent of a real number. 



DAP Calling 

Sequence 



CALL TANH 

DAC ARG1 (a real number) 

(Return) 



FORTRAN Reference TANH(R) 



Method 



TANH = (e**(2*X)-l)/(e**(2*X)+l), where X = ARG1. 



Data Type of 
Arguments and 
Results 



This tangent function with a real argument results in a real number. 



Other Routines L$22, EXP, A$22, H$22, D$22 

Used 
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SECTION V 
COMPILER SUPPORT SUBROUTINES 



This section describes the compiler support subroutines, i.e., those subroutines which 
are not normally explicitly called by the FORTRAN programmer. These subroutines per- 
form conversions between data types, logical relationals, arithmetic operations, and miscel- 
laneous functions. 
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Purpose 



To add or subtract real numbers. (This subroutine has a high- 
speed version, A$22X. ) 



DAP Calling 
Sequence 



CALL A$22 
DAC ARG2 
(Return) 



(or S$22) 

(a real number) 



Method 



A$22 (Add) - The contents of ARG2 are added to the contents of the 
A- and B- registers after both numbers are unpacked and scaled. 
The result is normalized and the characteristic is adjusted. 

S$22 (Subtract ) - The value contained in ARG2 is negated and the 
add routine, A$22, is entered. 



Data Type of 
Arguments and 
Results 



< implicit real argument> ± < real argument> -*- < real result> 



Error Messages 



The error message "SA" is reported if an arithmetic overflow 
occurs, i. e. , the result is > 2**127. An undefined result is 
returned. 



Other Routines 
Used 



ARG$, N$22, F$ER 
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Purpose 



To add or subtract real numbers. (This routine requires the High- 
Speed Arithmetic Option.) 



DAP Calling 
Sequence 



CALL A$22X 
DAC ARG2 
(Return) 



(A$22, S$22 or S$22X) 
(a real number) 



Method 



A$22 (Add ) - The contents of ARG2 are added to the contents of the 
A- and B- registers after both numbers are unpacked and scaled. 
The result is normalized and the characteristic is adjusted. 

S$22 (Subtract ) - The value contained in ARG2 is negated and the 
add routine, A$22, is entered. 



Data Type of 
Arguments and 
Results 



< implicit real argument> ± < real argument> -«- < real result> 



Error Messages 



The error message "SA" is reported if an arithmetic overflow 
occurs, i. e. , the result is > 2**127. An undefined result is 
returned. 



Other Routines 
Used 



N$22, F$ER 
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Purpos 



To add a real argument to a complex number. 



DAP Calling 
Sequence 



CALL A$52 

DAC ARG2 (a real number) 

(Return) 



Method 



The following is the algorithm used to compute the operation of 
adding a real argument (ARG2) to the contents of the complex 
accumulator (Y): 

Y + ARG2 = A + B * I + ARG2 

= (A + ARG2) + B * I) 
where Y = A + B * I 



Data Type of 
Arguments and 
Results 



< implicit complex argument> + < real argument> -»- < complex 
re suit > 



Other Routines 
Used 



F$AT, H$55, L$22, A$22, H$22, L$55 
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Purpose 



To add complex numbers. 



DAP Calling 
Sequence 



CALL A$55 

DAC ARG2 (a complex number) 

(Return) 



Method 



The following is the algorithm used in the addition of two complex 
numbers (the contents of ARG2 and the complex accumulator): 

X+ARG2 = (A+B*I) + (M+N*I) = (A+M) + (B+N) * I 
where X = A+B*I and ARG2 = M+N*I 



Data Type of 
Arguments and 
Results 



< implicit complex argument> + < complex argument> — *- 

< complex re suit > 



Other Routines 
Used 



F$AT, H$55, SUB$, L$22, A$22, H$22, L$55 
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Purpose 



To add a real number to a double-precision number. 



DAP Calling 
Sequence 



CALL A$62 

DAC ARG2 (a real number) 

(Return) 



Method 



This subroutine calls DBLE to convert the real argument to a 
double-precision number and calls A$66 to perform the double- 
precision addition. 



Data Type of 

Arguments and 
Results 



< implicit double -precision argument> + < real argument> -*- 

< double -precision result> 



Other Routines 
Used 



F$AT, H$66, DBLE, A$66 



5-6 



AM74 



Purpos 



To add, subtract, multiply, or divide normalized, double- 
precision numbers. (This subroutine has a high-speed version, 
A$66X. ) 



DAP Callin 



£ 



Sequence 



CALL A$66 
DAC ARG2 
(Return) 



(or S$66, M$66, or D$66) 
(a double-precision number) 



Method 



Data Type of 
Arguments and 
Results 



The contents of ARG2 are added to, subtracted from, multiplied 
by, or divided into the contents of the double -precision accumu- 
lator (X). 

Add (A$66) - The numbers are unpacked and scaled to coincident 
places. The addition process takes place (X+ARG2), and the result 
is normalized. 

Subtract (S$66) - The numbers are unpacked and scaled to coincident 
places. The subtraction process takes place (X-ARG2), and the 
result is normalized. 

Multiply (M$66) - X*ARG2 = (X*2**E1) * (Y*2**E2) 

= X*ARG2*2** (E1+E2) 
Let X = (A+B*2** (-N)) 
and ARG2 = (C+D*2** (-N)) 

X*ARG2 = A*C + ((A*D+B*C) * 2**(-N)) 

The term B*D*2**(-2N) is ignored. 

The least significant bits of the product are: 

L*(A*C)+H*(A*D)+H*(B*C) 

Divide (D$66) - The quotient X/ARG2 is obtained by the binomial 
expansion of 1/X = X**(=l). The high-order and low-order parts 
(H and L) of the quotient are computed as follows: 

(A+B*2**(-N))/(C+D*2**(-N)) = (A+B-A*D/C)/C 

H = (A+B-A*D/C)/C 
L = remainder (H)/C 



< implicit double -precision argumetit> 



,, < double -precision 



argument> — ■+- < double -precision result> 



Error Messages 



The error message "AD" is printed if an addition or 
subtraction over /underflow occurs. 

The error message "PZ" is printed if a division by zero 
is attempted. 

The error message "MD" is printed if a multiplication or 
division over /underflow, occurs. 
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A$66 cont. 



After an error message is reported, the double-precision accumu- 
lator is loaded with the maximum ((2**128)- 1) or minimum 
(2**(-128)) value (as determined by the correct sign) before 
returning to the calling program. 



Other Routines N$66, F$ER, H$66, L$66, ARG$, AC1, AC2, AC3 

Used 



AM74 



Purpose 



To add, subtract, multiply, or divide normalized, double-precision 
numbers. (This routine requires the High-Speed Arithmetic Option. ) 



DAP Calling 
Sequence 



CALL A$66X (or A$66, S$66, S$66X, M$66, M$66X, D$66, D$66X) 

DAC ARG2 (a double-precision number) 

(Return) 



Method 



The contents of ARG2 are added to, subtracted from, multiplied by, 
or divided into the contents of the double-precision accumulator. See 
A$66, described on the preceding pages, for a detailed description of 
the methods used. 



Data Type of 
Arguments and 
Results 



f+1 



< implicit double -precision argument > 



,. ■< double -precision 



argument> -*- < double -precision result> 



Error Messages See Error Messages for A$66. 



Other Routines N$66, F$ER, H$66, L$66, ARG$, AC1, AC2, AC3 



Used 
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A$81 



P urpose To add an integer value (I) to the characteristic of the variable in the 

double-precision accumulator (effectively, multiplication by 2 1 ). 

DAP Calling CALL A$81 

Sequence DAC ARG2 (an integer value) 

(Return) 

Met hod The characteristic (base 2) of the value in the double-precision accum- 

ulator is increased (or decreased) by an integral value, ARG2. For 
example, if ARG2= 2 and the value in the double -precision accumu- 
lator is 8.0 (2 3,0 ), the result of this call would be 2 ' or 2 ' = 
32.0 (8.0*2 2 ). If the absolute value of the result is less than 2**(-128), 
a value of zero is returned. 

Data Type of < implicit double -precision argument> * (2**< integer argument>)-» 

Arguments and < double -precision result> 

Results 



Error Messages If there is exponent overflow, an "EQ" error message is reported and 

external locations AC1 and AC2 are loaded with the maximum value 
possible' ((2**128)- 1) with the sign of ARG2. 

Other Routines N$22, F$ER, AC1, AC2 

Used 
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AC1 

(AC2, AC3, AC4, AC5) 



Purpose To assign locations to be used as a double-precision or complex 

accumulator by the FORTRAN library routines. 

U se AC1, AC2, AC3: double-precision accumulator. 

AC1, AC2: complex accumulator, real portion. 

AC3, AC4: complex accumulator, imaginary portion. 

AC5: error flag. 
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Purpose To convert the indirect address of an argument to its corresponding 

direct address. 



DAP Calling CALL ARG$ 

Sequence DAC* ARG2 (usually a subroutine entry) 

(Return) 

Method The address of the argument is returned in the index register. This 

subroutine may be used upon entering a subroutine to set up the return 
addre ss. 
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Purpose To convert an integer to a real number. 

DAP Calling CALL C$12 

Sequence (Return) 

Method Th e integer value in the A-register is placed in the B-register and 

the A-register is set to 045600 (octal), representing a characteristic 
such that the number fits the description given for a real number ex- 
cept that it is not "normalized. " A$22 (with argument = 
(040000, 000000), also unnormalized) is called to normalize the 
result. 



Data Type of The integer value in the A-register is converted to a real number and 

Arguments and placed in the A- and B-registers. 

Results 



Other Routines A$22, N$22 

Used 
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Purpose To convert an integer to a double-precision number. 

DAP Calling CALL C$16 

Sequence (Return) 

Method The integer in the A-register is normalized and converted to real by 

calling C$12. This real value is then converted to a double-precision 
number by calling C$26. The result is placed in the double-precision 
accumulator. AC 1 contains the contents of the B-register (the real 
exponent), AC2 contains the contents of the A-register (the most 
significant word of the fraction), and AC3 contains a word of zeros. 

Data Type of The integer value in the A-register is converted to a double-precision 

Arguments and number and placed in the double-precision accumulator. 

Results 

Other Routines C$12, C$26 

Used 



5-14 AM74 



Purpose 



To convert a real number to an integer. 



DAP Calling 
Sequence 



CALL C$21 
(Return) 



Method 



This subroutine scales the real number in the A- and B-registers to 
23 bits by adding the octal value 045700 (2**22) to truncate the 
fractional part of the real number. The result is in the A-register. 



Data Type of 
Arguments and 
Results 



The real number in the A- and B-registers is converted to an integer 
and returned in the A-register. 



Error Messages 



The message "RI" is reported if the integer (I) is too large when con- 
verted from real to integer. The integer must be in the following 
range: -2 15 < I < 2 15 -1. An undefined result is returned in the A- 
register. 



Other Routines 
Used 



N$22, A$22, F$ER 



5-15 



AM74 



Purpose 



To convert a real number to a complex number. 



DAP Calling 
Sequence 



CALL C$25 
(Return) 



Method 



The A- and B- registers are stored in AC1 and AC2, respectively 
(the real part of the complex number), and AC3 and AC4 (the imagin- 
ary part of the complex number) are set to zeros. 



Data Type of 
Arguments and 
Results 



The real argument in the A- and B- registers is converted to a com- 
plex number and stored in the complex accumulator (AC1, AC2, AC3, 
and AC4). 



Other Routines 
Used 



H$22, CMPLX 



5-16 



AM74 



Purpose 



To convert a real number to a double -precis ion number. 



DAP Calling 
Sequence 



CALL C$26 
(Return) 



Method 



The number in the A- and B-registers is placed in AC1 and AC2. 
AC3 is cleared and the routine exits. 



Data Type of 
Arguments and 
Results 



The real number in the A- and B-registers is converted to double- 
precision and placed in the double-precision accumulator. 



Other Routines 
Used 



AC1, AC2, AC3 



5-17 



AM74 



Purpose 



To convert a double-precision number to an integer. 



DAP Calling 
Sequence 



CALL C$61 
(Return) 



Method 



This subroutine calls C$62 to convert the number in the double - 
precision accumulator to real and calls C$21 to convert the real 
number to integer. 



Data Type of 
Arguments and 
Results 



The double-precision value in the double-precision accumulator is 
converted to an integer and placed in the A-register. 



Other Routines 

Used 



C$62, C$21 



5-1! 



AM74 



Purpose To convert a double-precision number to a real number. 

DAP Calling CALL C$62 or CALL SNGL 

Sequence (Return) DAC ARG1 (a double-precision 

(Return) number) 

Method AC1 and AC2 (the exponent and the most significant part of the 

fraction of the number in the double-precision accumulator) or the 
first two words of ARG1 (if SNGL is called) are loaded into the A- 
and B- registers. The least significant part of the fraction (AC3. 
word 3) is not considered in the result. 



Data Type of The double-precision value in the double-precision accumulator or 

Arguments and in ARG1 is converted to a real number and placed in the A- and B- 

Results registers. 

Other Routines L$22, N$66, N$22, L$66, AC1, AC2 
Used 



5-19 AM74 



Purpose 



To convert the exponent of the value in the double-precision accumu- 
lator to an integer. 



DAP Calling 
Sequence 



CALL C$81 
(Return) 



Method 



Extract the characteristic (base 2) from the value in the double- 
precision accumulator (AC1) and convert it to an integer. 



Data Type of 
Arguments and 
Results 



The characteristic of the double-precision argument is converted to 
an integer. 



Other Routines 
Used 



AC1 



5-20 



AM74 



D$11 



Purpose To divide two integers. (This subroutine has a high-speed version, 

~~~ D$UX. ) 

DAP Calling CALL D$ll 

Sequence DAC ARG2 (integer divisor) 

(Return) 



Method The numerator (an integer value) should be in the A-register upon 

entrance to this subroutine. If the denominator, ARG2, is zero, 
an overflow occurs and an error message is reported. If both 
arguments are nonzero, the numerator is positioned in the A- and 
B-registers and the division is performed. The results are ex- 
amined for the special case (-32, 768/-1) which is treated as an over- 
flow. If the results are in the range of -32, 768 to + 32, 767, D$ll 
returns to the calling program with the quotient in the A-register 
and the remainder in the B-register. The integer answer is in the 
A-register. 



Data Type of < implicit integer argument> / < integer argument> -+■ < integer 

Arguments and re suit > 

Results 

Error Messages The error message "IZ" is reported if a division by zero is at- , 

tempted. The maximum value is output (-32,768 if negative or 
+ 32, 767 if positive). A division of -32, 768 by -1 also causes "IZ" 
to be reported; D$ll returns a value of + 32, 767, the maximum 
value possible. 



Other Routines ARG$, F$ER 

Used 



5-21 AM74 



DS11X 



Purpose 



To divide two integers. (This routine requires the High-Speed 
Arithmetic Option. ) 



DAP Calling 
Sequence 



CALL D$11X (or D$ll) 

DAC ARG2 (integer divisor) 

(Return) 



Method 



See "Method" for D$ll. 



Data Type of 
Arguments and 
Results 



< implicit integer argument> / < integer argument> — •- < integer 
re suit > 



Error Messages 



See "Error Messages" for D$ll, 



Other Routines 
Used 



ARG$, F$ER 



5-22 



AM74 



Purpose To divide two real numbers. (This subroutine has a high- 

speed version, D$22X. ) 



See M$22. 



5 " 23 AM74 



Purpose 



To divide two real numbers. (This subroutine requires the High- 
Speed Arithmetic Option. ) 



DAP Calling 
Sequence 



CALL D$22X (or D$22) 

DAC ARG2 (the real divisor) 

(Return) 



Method 



Data Type of 
Arguments and 
Results 



This subroutine divides the real number in the A- and B-registers 
(X) by the real argument, ARG2 (Y). The division is performed by 
multiplying X by the reciprocal of Y, i.e., X*l/Y. Newton's 
method for 1 /Y is: 

R(l) = R(0) * (2-R(0)*Y 

where 

R(0) - 1/H(Y), H(Y) being the high-order 15 bits of Y 

X * (1/Y) = X * R(l) = X * R(0) * (2-R(0)*Y) 
< implicit real argument> / < real argument> -*■ < real result> 



Error Messages 



A "DZ" error message is typed if division by zero is attempted. 
A value of is returned if the dividend is also 0. The signed 
maximum value (±1. 7E38) is returned if the dividend is nonzero. 

An "SM" error message is reported if an arithmetic overflow 
occurs. The signed maximum value (±1.7E38) is returned. 

A value of is returned for an overflow. 



Other Routines 

Used 



N$22, F$ER 



5-24 



AM74 



Purpose To divide a complex number by a real number. 



DAP Calling CALL D$52 

Sequence DAC ARG2 (a real number) 

(Return) 



Method This subroutine divides the complex value in the complex 

accumulator (Y) by the real argument, ARG2. 

Y/ARG2 = (A+B*I)/ARG2 = A/ARG2+B*I, where Y = A+ B*I 

Data Type of < implicit complex argument> /< real argument> -*- <complex 

Arguments and re suit > 

Results 



Other Routines F$AT, H$55, SUB$, L$22, D$22, H$22, L$55 

Used 



5 " 25 AM74 



Purpose 



To divide two complex numbers. 



DAP Calling 
Sequence 



CALL D$55 

DAC ARG2 (complex divisor) 

(Return) 



Method 



The following algorithm is used to compute the operation of 
dividing two complex numbers. The contents of the complex 
accumulator (X) are divided by the contents of ARG2 (Y). 

X/Y = (A+B*I)/(M+N*I) 

where X - A+ B*I and Y = M+ N*I 

= (A+ B*I)*(M-N*I)/(M+N*I)*(M-N*I) 

= (A+ B*I)*(M-N*I)/(M**2+N**2) 

= (A*M+B*N+ B*M*I-A*N*I)/(M**2+N**2) 

= (A*M+ B*N)/(M**2+ N**2)+ (B*M*I-A*N*I)/(M**2+ N**2) 

= (A*M+ B*N)/(M**2+ N**2)+ (I*(B*M-A*N))/(M**2+ N**2) 



Data Type of 
Arguments and 
Results 



< implicit complex argument> / < complex argument> 
re suit > 



< complex 



Other Routines 
Used 



F$AT, H$55, SUB$, L$22, M$22, H$22, A$22, D$22, S$22, 
N$22, L$55 



5-26 



AM74 



Purpo se 



To divide a double-precision number by a real number. 



DAP Calling 
Sequence 



CALL D$62 
DAC ARG2 

(Return) 



(a double-precision number) 



Method 



This subroutine calls DBLE to convert the real divisor (ARG2) 
to a double-precision number and calls the double-precision 
divide routine (D$66). 



Data Type of 
Arguments and 
Results 



< implicit double -precision argument > / < real argument > 

< double -precision result> 



Other Routines F$AT, H$66, DBLE, L$66, D$66 



Used 



5-27 



AM74 



Purpose To divide normalized double -precision numbers. 



See A$66. 



5-2 8 AM74 



Purpose 



To calculate the value of an integer raised to an integer power. 
(This subroutine has a high-speed version, E$11X.) 



DAP Calling 
Sequence 



CALL E$ll 
DAC ARG2 

(Return) 



(the integer exponent) 



Method 



The implicit integer argument in the A-register and the integer 
exponent, ARG2, are first examined for the combinations listed 
below. If one of these combinations is found, the answer is loaded 
in the A-register for return to the calling program. 



Value in A -I 


Register 


Exponent 


Answer 


I 









1 












1 









- 


+32767 









+ 





1 






J 


1 


-1 






even 


1 


-1 






odd 


-1 


1 






- 






Otherwise, the value of the expression is calculated and returned 
in the A-register. The maximum or minimum value computed 
may not exceed +32, 767 or -32, 768. 



Data Type of 
Arguments and 
Results 



< implicit integer argument> ** < integer argument> — ►- < integer 
re suit > 



Error Messages 



The error message "II" is reported and +32,767 is returned if 
overflow occurs or if I = and J is negative (1/0). The value 
-32,768 is returned if I<-2, J is odd, and overflow occurs. 



Other Routines 
Used 



ARG$, M$ll, F$ER 



5-29 



AM74 



Purpose 



To calculate the value of an integer raised to an integer power. 
(This subroutine requires the High-Speed Arithmetic Option.) 



DAP Calling 
Sequence 



CALL E$11X (or E$ll) 

DAC J (the integer exponent) 

(Return) 



Method 



See "Method" for E$ll, 



Data Type of 
Arguments and 
Results 



< implicit integer argument> **< integer argument> -»- < integer 
result> 



Error Messages See "Error Messages" for E$l 1. 



Other Routines 
Used 



ARG$, F$ER 



5-30 



AM74 



Pu 



rpose 



To calculate the value of a real number raised to an integer power. 



DAP Calling 
Sequence 



CALL E$21 

DAC ARG2 (the integer exponent) 

(Return) 



Method 



A**ARG2 is evaluated by multiplying A by itself ARG2-1 times. 
The sign is determined by the sign of the number in the A- and 
B- registers and whether I is odd or even. 



Data Type of 
Arguments and 
Results 



< implicit real argument> **< integer argument> -*- < real result> 



Other Routines 
Used 



ARG$, M$22, D$22 



5-31 



AM74 



Purpose 



To calculate the value of a real argument raised to a real power. 



DAP Calling 
Sequence 



CALL E$22 
DAC ARG2 

(Return) 



(the real exponent) 



Method 



X**ARG2 is evaluated as e**(ARG2*log(X)). 



Data Type of 
Arguments and 
Re suits 



< implicit real argument> **< real argument> -*■ < real result> 



Other Routines 
Used 



ARGS, ALOG, M$22, EXP 



5-32 



AM74 



Purpose To calculate the value of a real number raised to a double -precision 

power. 



DAP Calling CALL E$26 

Sequence DAC ARG2 (the double-precision exponent) 

(Return) 



Method B**ARG2 is evaluated as e**ARG2*log(B)). 



Data Type of < implicit real argument> ** < double -precision argument> 

Arguments and < double-precision result> 

Results 



Other Routines F$AT, C$26, H$66, DLOG, M$66, DEXP 

Used 



5-33 AM74 



Purpose To calculate the value of a complex quantity raised to an integer 

power. 

DAP Calling CALL E$51 

Sequence DAC ARG1 (the integer exponent) 

(Return) 



Method The number in the complex accumulator is multiplied by itself 

.. _ ARG1-1 times. 



Data Type of <implicit complex argument> ** < integer argument> -»- < complex 

Arguments and re suit > 

Results 



Other Routines F$AT, H$55, IABS, L$55, M$55, D$55 

Used 



5-34 AM74 



Purpose To calculate the value of a double-precision number raised to an 

integer power. 

DAP Calling CALL E$61 

Sequence DAC ARG2 (the integer exponent) 

(Return) 

Method This routine checks for an even -numbered exponent, squares the 

number in the double-precision accumulator, and divides the integer 
argument (the exponent) by 2 until the exponent divided by 2 = 1. 
If the exponent is odd, the computed value (Dl""*) is multiplied by 
the original double-precision number before exiting. 

Data Type of < implicit double -precision argument> ** < integer argument> — »- 

Arguments and < double -precision result> 

Results 



Other Routines F$AT, H$66, L$66, D$66, D$ll, M$ll, M$66 

Used 



5-35 AM74 



Pur pose To calculate the value of the number in the double-precision accumu- 

lator raised to a real power. 

DAP Calling CALL E$62 

Sequence DAC ARG2 (the real exponent) 

(Return) 

Method B**ARG2 is evaluated as e**(ARG2*DLOG(B)), where B = the 

contents of the double -precision accumulator. 

Data Type of < implicit double -precision argument> ** < real argument> -*- 

Arguments and < double -precision result> 

Results 

Other Routines F$AT, H$66, DLOG, M$62, DEXP 

Used 



5-36 AM74 



Pu rpo s e 



To calculate the value of a double-precision value raised to a 
double-precision result. 



DAP Calling 
Sequence 



CALL E$66 
DAC ARG2 

(Return) 



(the double-precision exponent) 



Method 



B**ARG2 is evaluated as e**(ARG2*LOG(B)), where B = the 
contents of the double-precision accumulator. 



Data Type of 
Arguments and 
Results 



< implicit double -precision argument > ** < double -precision argument > 
-»- < double -precision result> 



Other Routines 
Used 



F$AT, H$66, DLOG, M$66, DEXP 



5-37 



AM74 



Purpose 



To store (hold) the contents of the A- and B -registers in memory. 



DAP Calling 
Sequence 



CALL H$22 
DAC ARG1 

(Return) 



(location in which the contents of the A- and B- 
registers are to be stored) 



Method 



The contents of memory at the location specified by the argument 
address, ARG1, are replaced by the contents of the A- and B- 
registers. The contents of the A- and B -registers remain un- 
changed. 



Data Type of 
Arguments and 
Results 



This subroutine stores a real number in the argument address. 



Other Routines 
Used 



ARG$ 



5-1,1 



AM74 



Purpose 



To hold (store) the contents of the complex accumulator in memory. 



DAP Calling 
Sequence 



CALL H$55 
DAC ARG1 

(Return) 



(location in which the contents of the complex 
accumulator are to be stored) 



Method 



The contents of memory at the location specified by the argument 
address, ARG1, are replaced by the contents of the complex accumu- 
lator. The contents of the accumulator remain unchanged. 



Data Type of 
Arguments and 
Results 



This subroutine stores a complex number in the argument address. 



Other Routines 
Used 



ARG$, AC1, AC2, AC3, AC4 



5-39 



AM74 



Purpose 



To hold (store) the contents of the double -precision accumulator 
in memory. 



DAP Calling 
Sequence 



CALL H$66 
DAC ARG1 

(Return) 



(location in which the contents of the double- 
precision accumulator are to be stored) 



Method 



The contents of memory specified by the argument address, ARG1, 
are replaced by the contents of the double-precision accumulator. 
The contents of the accumulator are unchanged. 



Data Types of 
Arguments and 
Results 



This subroutine stores a double -precision number in the argument 
address. 



Other Routines 
Used 



ARG$, AC1, AC2, AC3 



5-40 



AM74 



Purpose 



To load a real number into the A- and B- registers. 



DAP Calling 
Sequence 



CALL L$22 or CALL REAL 
DAC ARG1 (a real number) 

(Return) 



Method 



This subroutine calls ARG$ to place the address of the argument, 
ARG1, into the index register. ARG1 is then loaded into the A- 
and B- registers. 



Other Routines 
Used 



ARG$ 



5-41 



AM74 



Purpose 



To form an INCLUSIVE OR from memory with the value in the 
A- register. 



DAP Calling 
Sequence 



CALL L$33 

DAC ARG1 (an integer value) 

(Return) 



Method 



The value in the A- register is EXCLUSIVELY ORed, ANDed, ant 
EXCLUSIVELY ORed again with the argument, ARG1. 



5-42 



AM74 



Purpose 



To load a complex number into the complex accumulatoi 



DAP Calling 
Sequence 



CALL L$55 

DAC ARG1 (a complex number) 

(Return) 



Method 



This subroutine calls ARG$ to place the address of the argument, 
ARG1, into the index register. ARG1 is then loaded into the comple 
accumulator. 



Other Routines ARG$, AC1, AC2, AC3, AC4 



Used 



5-43 



AM74 



Purpose 



To load a double-precision number into the double-precision 
accumulator. 



DAP Calling 
Sequence 



CALL L$66 

DAC ARG1 (a double-precision number) 

(Return) 



Method 



This subroutine calls ARG$ to place the address of the argument, 
ARG1, into the index register. ARG1 is then loaded into the 
double-precision accumulator. 



Other Routines 
Used 



ARG$, AC1, AC2, AC3 



5-44 



AM74 



M$11 



Purpose 



To multiply two integers, 
M$11X.) 



(This subroutine has a high-speed version, 



DAP Callin 



g 



Sequence 



CALL M$ll 
DAC ARG2 

(Return) 



(integer multiplier) 



Method 



This subroutine multiplies the value in the A- register by the integer 
argument, ARG2. If either or both are negative, a sign counter is 
incremented and the negative value(s) are made positive. The multi- 
plier, ARG2, is loaded into the B- register and shifted to place the 
low- order bit of the multiplier in the C- register. The C-bit is 
tested and if it is set, the multiplicand is added to the A- register. 
The A- and B- registers are shifted together 1 bit, with the new low- 
order bit going into the C- register, and so forth, for 16 shifts. 
When these right shifts are completed, the bits are shifted back into 
the A- register, one at a time, checking for overflow. The positive 
or negative result is returned in the A- register. 



Data Type of 
Arguments and 
Results 



< implicit integer argument> *< integer argument> -*- < integer result> 



Error Messages 



When an over /underflow occurs, the error message "IM" is reported. 
The subroutine returns with +32, 767 in the A- register if the answer 
is positive, or -32,768 if it is negative. 



Other Routines 
Used 



ARG$, F$ER 



5-45 



AM74 



M$11X 



Purpose 



To multiply two integers. (This subroutine requires the High-Speed 
Arithmetic Option. ) 



DAP Calling 
Sequence 



CALL M$11X (or M$ll) 

DAC ARG2 (an integer value) 

(Return) 



Method 



This subroutine multiplies the value in the A- register by ARG2. 
The result is then examined for over/ underflow (see "Error 
Messages"). If the result is in the proper range, the signed result 
is returned to the calling program in the A- register. 



Data Type of 
Arguments and 
Results 



<implicit integer argument > * <integer argument > — »- <integer result > 



Error Messages 



See "Error Messages" for M$ll. 



Other Routines 
Used 



ARG$, F$ER 



5-46 
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Purpose 



To multiply or divide two real numbers. (This subroutine has a 
high-speed version, M$22X. ) 



DAP Calling 
Sequence 



CALL M$22 (or D$22) 
DAC ARG2 (multiplier 

(Return) or divisor) 



The multiplicand (M$22) or dividend 
(D$22) must be in the A- and B- 
registers. The sign, exponent, and 
most significant bits will be in the 
B- register. 



Method 



X*Y = (X*2**B)*(Y*2**C), where X = the value in the A- and B. 

registers 

Y = ARG2 

= ABS(X)*ABS(Y)*2**(B+C) 
ABS(X)*ABS(Y) = X(1)*Y(1)+(X(1)*Y(2)+X(2)*Y(1))*2**-15 
The most significant part of the product is H(X(1)*Y(1)) and 
the least significant part is L(X(1)* Y(1))+H(H(1)*Y(2))+H(X(2)* 
Y(l))*2**-15. 

Newton's method for 1/Y is R(l) = R(0)*(2-R(0)# Y), where 
R(0) = 1/H(Y), H(Y) being the high-order 15 bits of Y. 
X(l/Y) . X*R(1) = X*R(0)*(2-R(0)*Y). 



Data Type of 
Arguments and 
Results 



< implicit real argument> *< real argument> -*- < real result> 



Error Mess 



ages 



Multiplication - If there is underflow, a value of zero is returned 
with no error message. 

If there is overflow, an "SM" error message is reported and the 
maximum value ((2**128)- 1) is returned in the A- and B-registers. 

Division - If division by zero is attempted, a "DZ" error message 
is reported and the result in the A- and B-registers is undefined. 

If the divisor is unnormalized, an "SD" error message is reported 
and the result in the A- and B-registers is undefined. 



Other Routines 
Used 



N$22, ARG$, F$ER 



5-47 



AM74 



Purpose 



To multiply two real numbers. 



DAP Calling 
Sequence 



CALL M$22 
DAC ARG2 

(Return) 



(a real number) 



Method 



Dat a Type of 

Arguments and 
Results 



X*Y 



(X*2**B)*(Y*2**C), where X = the value in the A- 

registers 
Y = ARG2 



and B- 



= ABS(X)*ABS(Y)*2**(BC) 
ABS(X)*ABS(Y) = X(1)*Y(1)*(X(1)*Y(2)+X(2)*Y{1))*2**-15 
The most significant part of the product is H((X( 1)*Y( 1)) and the 
least significant part is L(X(1)* Y{ 1))+H(H( 1)*Y(2))+H(X(2)* Y( 1)) 

*2**-15, where H(X(1)*Y(1)) is the most significant part of the 
product X(1)*Y(1) and L(X(1)*Y(1)) is the least significant part 
of that product. 
< implicit real argument> * < real argument> -»- < real result> 



Error Messages 



Underflow - A value of zero is returned with no error message. 
Overflow - An "SM" error message is reported and a signed maxi- 
mum value (±1.7E38) is returned. 



Ot her Routines 
Used 



F$ER 



5-4£ 
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Purpose 



To multiply a complex number by a real number. 



DAP Calling 
Sequence 



CALL M$52 

DAC ARG2 (a real number) 

(Return) 



Method 



Y*X = (A+B*I)*X = A*X+(B*X)*I 
where Y = A+B*I (in the complex accumulator) 
X = ARG2 



Data Type of 
Arguments and 
Results 



< implicit complex argument> *< real argument> -*- < complex result> 



Other Routines F$AT, H$55, SUB$, L$22, M$22, H$22, L$55 

Used 



5-49 
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Purpose 



To multiply complex numbers. 



DAP Calling 

Sequence 



CALL M$55 

DAC ARG2 (a complex value) 

(Return) 



Method 



This routine multiplies the contents of the complex accumulator (X) 

by the value in ARG2 (Y). 

X*Y = (A+B*I) (M+N*I) = A*M-B*N+(A*N+B*M)*I 

where X = A+B*I and Y = M+N*I. 



Data Type of 
Arguments and 
Results 



< implicit complex argument> * < complex argument> -*- < complex 

re suit > 



Other Routines 
Used 



F$AT, H$55, SUB$, L$22, M$22, H$22, S$22, N$22, A$22, L$55 
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Purpose 



To multiply a double-precision number by a real number. 



DAP Calling 
Sequence 



CALL M$62 

DAC ARG2 (real multiplier) 

(Return) 



Method 



This subroutine calls DBLE to convert the real multiplier to a 
double-precision number and calls the double-precision multiply 
routine (M$66). 



Data Type of 
Arguments and 
Results 



< implicit double -precision argument> * < real argument> -*- 

< double -precision result> 



Other Routines F$AT, H$66, DBLE, M$66 

Used 
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Purpose To multiply normalized, double-precision numbers. 

See A$66. 
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Pur P ose To determine the TWOs complement of a real number. 

DAP Calling CALL N$22 

Se( l uence DAC ARG1 (a real number) 

(Return) 

Meth ° d The C " bit is preset on entrance to this routine to provide a true TWOs 

complement if the low- order word is found to be zero. The C-bit is 
reset when this is not the case, and the A- and B- registers are TWOs 
complemented normally. 

Data Type of The TWOs complement of the real argument is computed and the 

Arguments and routine exits with the the real result in the A- and B- registers 

Results 6 
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Purpos 



To obtain the complement of a logical value. 



DAP Calling CALL N$33 

Sequence (Return) 

Method The least significant bit of the argument in the A- register is 

logically complemented, changing its value from true to false 

(1 to 0) or false to true (0 to 1). 
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Purpose 



To negate a complex quantity. 



DAP Calli 



mg 



Sequence 



CALL N$55 
(Return) 



Method 



The signs of the real part and the complex part of the complex 
number are negated. The result is in the complex accumulator. 



Data Type of 
Arguments and 
Results 



The complex argument is negated and the subroutine exits, with the 
complex result in the complex accumulator. 



Other Routines H $55, SUB$, L$22, N$22, H$22, L$55 



Used 
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Purpose 



To negate a double-precision number. 



DAP Calling 
Sequence 



CALL N$66 
(Return) 



Method 



This subroutine negates the value in the double-precision accumu- 
lator. The double-precision word is effectively TWOs complemented 
as follows: 



1. 



The lowest order word, AC3, word 3, is tested for zero. If it 
is not zero, the word is TWOs complemented. If it is zero, the 
C-bit is set. 

2 AC2, word 2, is tested for zero. If it is not zero, the word is 
ONEs complemented and the C-bit is added. If it is zero and the 
C-bit is set, no action is taken. If the C-bit is not set, the word 
is ONEs complemented. 

3. AC1, word 1, is ONEs complemented, and the C-bit, if set, is 
added. The negated result is left in the double-precision accumu 
lator. 



Data Type of 
Arguments and 
Results 



The double-precision argument is negated and the routine exits with i 
a double-precision result in AC 1, AC2, and AC 3. 



Other Routines 
Used 



AC1, AC2, AC3 
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Purp ose To subtract real numbers. (This subroutine has 

a high-speed version, S$22X. ) 

See A$22. 
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Purpose To subtract real numbers. (This subroutine requires the 

High-Speed Arithmetic Option. ) 

See A$22X. 
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Purpose 



To subtract a real number from a complex number to obtain a com- 
plex result. 



DAP Calling 
Sequence 



CALL S$52 

DAC ARG2 

(Return) 



(a real number) 



Method 



Y-X = A+B*I-X = (A-X)+B*I 
where Y = A+B*I, X = ARG2 



Data Type of 
Arguments and 
Results 



< implicit complex argument> -<real argument> -*- < complex result: 



Other Routines F$AT, H$55, L$22, S$22, H$22, L$55 



Used 



5-59 



AM74 



Purpose 



To subtract two complex numbers. 



DAP Calling 
Sequence 



CALL S$55 

DAC ARG2 (the complex subtrahend) 

(Return) 



Method 



This subroutine subtracts ARG2(Y) from the value in the complex 

accumulator (X): 

X - Y = (A+B*I) - (M+N*I) = A*M+(B*N)*I 

where X = A+B*I, Y = M+N*I 



Data Type of 
Arguments and 
Results 



<implicit complex argument> -<complex argument> 
re suit > 



< complex 



Other Routines 
Used 



F$AT, H$55, SUB$, L$22, S$22, N$22, H$22, L$55 
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Purpose 



To subtract a real argument from a double -precision number. 



DAP Calling 
Sequence 



CALL S$62 

DAC ARG2 

(Return) 



(a real number) 



Method 



This subroutine calls DBLE to convert the real argument to a double - 
precision number and enters the double -precision subtraction 
routine (S$66). 



Data Type of 
Arguments and 
Results 



< implicit double -precision argument> — < real argument> -+■ 

< double -precision result> 



Other Routines 
Used 



F$AT, H$66, DBLE, S$66, N$66 
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Purpose To subtract normalized, double-precision numbers. 

See A$66. 
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SNGL 



Purpose To convert a double -precis ion number to a real number. 

See C$62. 
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SUBS 



Purpose To calculate the address of a referenced array element or to calcu- 

late the array size. 

DAP Calling CALL SUB$ 

Sequence DAC or DAC* ARRAY TABLE1 

DAC or DAC* SUBSCRIPT 1 

DAC or DAC* SUBSCRIPT 2 



DAC or DAC* SUBSCRIPT N 
(Return) 



CALL SIZ$ 

DAC or DAC* ARRAY TABLE2 

(Return) 



Method ARRAY TABLE1 Layout 

DAC or DAC* ARRAY 

OCT L (number of words per array element) 

DEC DIMENSION 1 

DEC DIMENSION 2 



DEC DIMENSION N 

OCT O (end of dimension list) 

ARRAY TABLE2 Layout 

DAC or DAC* ARRAY 
OCT KEY 

OCT or DAC* DIMENSION 1 
OCT or DAC* DIMENSION 2 



OCT or DAC* DIMENSION N or OCT ARRAY SIZE 

or omitted 

The KEY bit pattern is CVDDDDDDDDDDDLLL, where 

C = - no array bounds checking 

C = 1 - array bounds checking 

V = - last word of array table is array size 

V = 1 - last word of array table is dimension 

If C = and V = 0, the last dimension word of the array table is 
omitted. 

D = dimensionality - limited to 2047 

L = number of words per array element 
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SUBS cont. 

Note that L is determined by the data type of the array as follows: 

Data Type of Array 

L = 1 - integer or logical 

2 - real 

3 - double-precision 

4 - complex 

Let S denote the array starting address, L the number of words per 
array element, S(I) the Ith subscript value, and D(I) the Ith dimension 
for an N- dimensional array A where N > 1. 

The address of the array element A(S(1), S(2), . . . S(N)) is given by 
S+ L(..., (S(N)-1)*D(N-1) +...+ (S(2)-1)*D(1) + (S(l)-D) 

Error Messages The error message "AO" (array overflow) is reported if the array 

element referenced is outside the bounds of the array. Only the final 
array element referenced is checked for legality, not individual sub- 
script values. 

Other Routines M$ll, F$ER 

Used 



5-65 AM74 



Z$80 



Purpose 



To clear (zero- out) the exponent of the variable in the double- 
precision accumulator. 



DAP Calling 
Sequence 



CALL Z$80 
(Return) 



Method 



Extract the value in AC1 and replace the characteristic (base 2) in 
bits 2-9 with zeros. 



Other Routines 
Used 



AC1 
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SECTION VI 
RUN -TIME AND CONTROL SUBROUTINES 



This section describes the routines which: control input and output by selecting and 
activating the proper device drivers; provide buffers; and edit and trace all I/O. 
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Purpose 



To transfer an array from or to the input or output data list. 



DAP Calling 
Sequence 



CALL 


F$AR 


OCT 


m 


DAC 


a 


(Return) 




or 




CALL 


F$Lx 


OCT 


m 


DAC 


a 


(Return) 





For use with DAP 
Number of words in array- 
Location of first word 



As used by the compiler 
x = 1 for integer 

= 2 for real 

= 3 for logical 

= 5 for complex 

= 6 for double precision 



FORTRAN 
Reference 



DIMENSION 1(5) 



WRITE (x, f) list 



ITOnTD ATM ctitcTYipnt nnr>-iV, «»■»•> 



Method 



Whenever the data list requires data from the internal source 
or data to be stored in the internal source, exit is made from 
F$10 to the next location in the data list. Elements of the data 
list may be variables, subscripted variables, or array names. 
Mode may be integer, real, logical, complex, or double 
precision. For easier data transmission, all list elements 
are assumed to be arrays; the mode is determined by 
appropriate format descriptors. For each item in the list, 
the three -word calling sequence (above) is generated. 



Other Routines 
Used 



F$10, F$CB, F$ER 
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Purpose 



To transfer a variable number of arguments from the 
calling routine to the called routine. 



DAP Calling 
Sequence 



DAC 
CALL 
DEC 
DAC 



F$AT 

n 

ARG1 



Entry point 

Number of arguments to be transferred 
Address at which first argument is stored 



DAC 
(Return) 



ARGn Address at which last argument is stored 



FORTRAN 
References 



CALL SUB1 (ARG1, ARG2, ..., ARGn) 

Where SUB1 is any subroutine and ARG1 
through ARGn are any constants, variables, 
arrays, etc. 



Method 



When arguments are to be transferred from a calling routine 
to a subroutine, a call to F$AT is generated by the compiler. 
The number of arguments specified by the first pseudo- 
operation following the call are transferred from the calling 
routine to the subroutine. All levels of indirect addressing 
are removed before an argument is transferred. ARG1 is 
the beginning location of the block into which the arguments 
are to be placed. 



Data Type 
of Arguments 



Arguments are direct relative addresses. 
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Purpose 



To connect the calling program with the magnetic tape 
rewind routine. 



DAP Calling 
Sequence 



CALL F$Bx x = 5, 6, 7, 8, 9 or the previously defined 

variable n (n = 5,6,7,8, or 9) 



FORTRAN 
Reference 



REWIND x 



x = 5, 6, 7, 8, 9 or the variable n 



Method 



This routine converts the logical magnetic tape unit number 
(5, 6, 7, 8 or 9) to the corresponding physical magnetic tape 
unit number (1,2, 3, 4, or 5) and then calls the REWIND 
routine to rewind the tape on that unit to the beginning of tape. 



Uther Routines 
Used 



C $MR 
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Purpose 



To close the buffers used for input or output. 



DAP Calling 
Sequence 



CALL 



F$CB 



FORTRAN 
Reference 



READ 

READ 

READ 

READ 

WRITE 

WRITE 

WRITE 



(n, m) list 

(n, m) 

(n) list 

(n) 

(n, m) list 

(n, m) 

(n) list 



n = device number; m 
statement number 



format 



Method 



At the end of the data list from a READ or WRITE statement, 
or when the format statement is exhausted (non-list), a call 
is issued to F$CB to close the buffer. The address of the 
buffer was determined by F$IO. F$CB checks for I/O mode 
and immediately closes the buffer if mode is input. If the 
buffer is formatted output, F$CB fills the remainder of the 
buffer with up to 134 spaces. If output is binary, the end of 
the buffer is filled with up to 120 zeros. 



Other Routines 
Used 



F$IO 
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Purpose 



To control the writing of an end-of-file mark on magnetic 
tape and a STOP code as an end-of-file on paper tape punch. 



DAP Calling 
Sequence 



CALL 


F$D2 


CALL 


F$Dx 


CALL 


F$Dn 



2 = paper tape punch 

x = logical tape unit number 5 through 9 

n = dummy device number and the A- 

register contains the value 2, 5, 6, 

7, 8, or 9 



FORTRAN 
References 



END FILE x 



x = 2, 5, 6, 7, 8, 9 or the variable n 



Method 



This routine converts the logical magnetic tape unit number 
(5, 6, 7, 8, or 9) to the corresponding physical magnetic tape 
unit number (1,2,3,4, or 5). It then calls the driver to 
write sn end-of-file mark on the specified magnetic tape. 
If x is 2, it calls the driver to punch a STOP code on the 
high-speed paper tape punch. 



Other Routines 
Used 



0$ME, 0$PS 
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Purpose 



To cause a mnemonic error indicator to be typed on the ASR-33 
when an object-time error is encountered in a specified routine. 



DAP Calling 
Sequence 



CALL F$ER ARG1 is the address of the indicator to 

DAC ARG1 be typed. The routine types the error 

indicator and halts if Sense Switch 3 
is not set. Pressing START after the 
halt causes the program to continue. 

If Sense Switch 3 is set, F$ER exits 
with no typeout and no halt. 



FORTRAN 
References 



CALL F$ER (2Hxx) xx = two ASCII characters to be typed 



Method 



The location of the object-time error mnemonic indicator 
is extracted from its relocatable address, AC5. AC5 is 
then loaded with the error mnemonic indicator and Sense 
Switch 3 is tested. If the switch is not set, a carriage 
return and line feed are issued. The error mnemonic is 
then printed, and the routine halts. If START is pressed 
at this point, a normal return is made to the calling routine. 
If Sense Switch 3 is set, return is made immediately to 
calling routine. 



Data Type 
of Argument 



The argument is the address of any two ASCII characters. 



Other Routines 
Used 



AC5, F$HT 
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Purpose 



To control back spacing of a record on magnetic tape. 



DAP Calling 
Sequence 



CALL F$Fx 



x = 5,6, 7, 8, 9 or the previously defined 
variable n. 



FORTRAN 
References 



BACKSPACE x x - 5, 6, 7, 8, 9 or the variable n. 



Method 



This routine converts the logical magnetic tape unit number 
(5, 6, 7, 8, or 9) to the corresponding physical magnetic tape 
unit number (1,2, 3,4, or 5) and then calls the driver to back 
space one record on the specified magnetic tape. 



Error 

Message 



An error message BF is reported if an end-of-file is 



encountered. 



Other Routines 
Used 



C$BR, F$ER 
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Purpose 



To process FORTRAN run-time assigned GO TO statements 



DAP Calling 
Sequence 



LDA 


PTR 


CALL 


F$GA 


DEC 


n 


DAC 


SI 


DAC 


S2 



Transfer address in A-register 

Number of statements in list 
Address of first statement 
Address of second statement 



DAC 



Sn 



Address of last statement 



FORTRAN 
Reference 



ASSIGN J TO I 
GO TO I, (K1.K2, 



J = statement number 
. Kn) I = integer variable name 

Ks = statement numbers 



Method 



This routine checks the address passed in the A-register against 
the statement address list that follows the call. If the address 
is found in that list, control passes to that statement. If not, 
a GO error is reported. No recovery from this error is 
possible. If the statement address list is empty (n = 0), the 
address is not checked. 



Data Type 
of Argument 



An address is passed to this routine in the A-register. 
trol is passed to the statement number at that address. 



Con- 



Other Routines 
Used 



F$ER 
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Purpose 



To process FORTRAN run-time computed GO TO statements. 



DAP Calling 
Sequence 



LDA 


PTR 


CALL 


F$GC 


DEC 


n 


DAC 


SI 


DAC 


S2 



Index to statement list 

Number of statements in list 
Address of first statement 
Address of second statement 



DAC 



Sn 



Address of last statement in list 



FORTRAN 
Reference 



GO TO (K1.K2, 



Kn),I I = integer value in the range 1 to n 

Ks = statement numbers 



Method 



The integer variable I or the content of PTR is treated as an 
index number; F$GC uses it to select the statement number 
from that position in the calling sequence. For example, 
if I = 3 then control is shifted to statement K3 or S3 in the 
above examples. If the index (I) is < 1 or >n, the computed 
GO TO statement is treated as a NOP. 



Data Type 
of Argument 



An integer value is passed to this routine in the A-register. 
Control is passed to the statement at the computed address. 
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FSHT 



Purpose 



To cause the computer to stop and print PA if a PAUSE 
statement has been encountered, or to print ST if a STOP 
statement has been encountered. 



DAP Calling 
Sequence 



CALL F$HT 

DAC '151724 Octal notation for ST 

CALL F$HT 

DAC '150301 Octal notation for PA 



FORTRAN 
Reference 



STOP or 

PAUSE 



Method 



This calling sequence is generated by the compiler when the 
STOP or PAUSE verb is encountered. The mnemonic ST or 
PA is placed in relocatable address AC5 and printed on the 
ASR. The A-register is restored and the program halts. 
Return to the calling program may be made by pressing START. 



Data Type 
of Argument 



The binary equivalent to the specified ASCII characters. 



Other Routines 
Used 



AC5 
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FSIO 



Purpose 



To perform input/output conversion, to edit input/output 
information, to accommodate the appropriate input /output 
device, and to provide buffers . 



DAP Calling 
Sequence 



DAC 


a 


CALL 


F$IO 


DAC 


BUF 


DAC* 


a 


CALL 


F$IO 


DAC 


BUF 



a = location of the format list from the 

READ statement 
BUF = buffer location 

The calling sequence for a FORTRAN 
WRITE with a and BUF as above 



FORTRAN 
References 



READ (n, x) list 
WRITE (n, x) list 
x FORMAT ( ) 



n = device number 

x = format statement number 



Method 



A FORTRAN READ /WRITE statement starts with a device 
number and a reference to a format statement, followed 
by an optional argument list. The first instruction generated 
by the READ/ WRITE statement is a coupling to the appro- 
priate device driver. The device driver then calls on 
F$IO passing the location of the format list, setting the entry 
location for the device driver and setting a flag to indicate 
input or output. F$IO then interprets the format list, 
character -by-character, taking whatever actions are 
required. Whenever data is required from or is to be stored 
in the internal source, exit is made from F$IO to the next 
location in the data list. 



Other Routines 
Used 



F$AR, F$CB, F$ER 
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F$m 



Purpose 



To control the typewriter keyboard input routine. 



DAP Calling 
Sequence 



CALL 

DAC n 

(Return) 



F$R1 



n - location of the format descriptor list 



FORTRAN 
References 



READ (l,f) list 
READ(l.f) 



f = FORTRAN statement number 



Method 



This subroutine connects the calling program with the 
I/O control subroutine (F$IO). Included in F$R1 is the driving 
logic needed to transmit input from the typewriter keyboard. 
When F$IO is called, the location of the format descriptor 
list (if any), the entry location of the driver subroutine, and 
a flag indicating input are transmitted. 

Whenever the F$IO subroutine requires data, return is made 
to the driver input entry of this subroutine, at which time up 
to 120 characters (terminated by a carriage return) are 
entered into the input buffer. 



Data Type 



Information is in ASCII format. 



Other Routines 
Used 



F$IO, I$AA 
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Purpose 



To control the paper tape input routine. 



DAP Calling 
Sequence 



FORTRAN 
References 



CALL 




F$R2 


DAC 




n 


(Return) 






READ (2, 


f) 


list 


READ (2, 


f) 




READ (2) 


list 


READ (2) 







Location of the format descriptor list 

(or 00000 if input is in binary format) 



Formatted READ where f is a format 
statement number 
Unformatted paper tape read 



Method 



This subroutine connects the calling program with the I/O 
control subroutine (F$IO). Included in F$R2 is the driving 
logic needed to transmit input from the paper tape reader. 

When F$IO is called the location of the format descriptor 

i;„i l;e -,„„\ <-!, „ „„t-^- r l^>r--a+-i^vi of +Vi e. Jrlvor siiKrniltitlp and 

a flag indicating input are transmitted. 

Whenever the F$IO subroutine requires data, it calls on the 
F$R2 driver subroutine, which assembles data into a 60-word 
buffer, three characters per word, if the input mode is 
binary, ii the input moue is L>^u, w uuaiciCi.^ j. o u.re 
assembled, two per word, into a 40 -word buffer. A carriage 
return character is replaced with as many blanks as needed 
to fill the rest of the input buffer. A tab character is 
replaced with as many blanks as needed to reach the next 
^redetermined tab position. 

A STOP code read in either binary or BCD mode causes the 
characters ST to be typed followed by a halt. Press START 
to continue. 



Data Type 



Information is in ASCII if formatted, or in binary if 
unformatted. 



Other Routines 
Used 



F$IO, I$PA, I$PB 
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Purpose 



To control the card input routine. 



DAP Calling 
Sequence 



FORTRAN 
References 



CALL 
DAC 




F$R3 

n 


(Return) 






READ (3, 
READ (3, 
READ (3) 
READ (3) 


f) list 

f) 
list 



Location of format descriptor list 
(00000 if input is binary) 



f = FORTRAN Statement number 
Unformatted read from a punch card 



Method 



This subroutine connects the calling program to the I/O 
control subroutine F$IO. Included in this subroutine is the 
driving logic required to input from the card reader. When 
F$IO is called, this subroutine transmits the location of the 
driver subroutine and a flag indicating input. 

When the F$IO subroutine requires data, it calls upon the 
F$R3 driver subroutine, which assembles 80 characters 
(two per word) into a 40 -word buffer, if mode is Hollerith 
(formatted), or into a 60 -word buffer in column binary- 
format (three words for every four columns) if binary. 
Return is made to F$IO where the buffer is processed. 



Data Type 



Information is in Hollerith if formatted, or in column 
binary if unformatted. 



Other Routines 
Used 



F$IO, I$CA, I$CB 
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F$R5-9 



Purpose 



To control reading of magnetic tape. 



DAP Calling 
Sequence 



CALL F$Rx 

DAC n 

(Return) 



x = 5, 6, 7, 8, or 9 

Location of the format descriptor list 

if formatted, or zero if unformatted 



FORTRAN 
References 



READ(x.f) list 

READ(x.f) 
READ(x) list 
READ(x) 



f = FORTRAN statement number and 

x = 5,6, 7, 8, or 9 
Unformatted read where x = 5, 6, 7, 8, or 9 



Method 



This subroutine connects the calling program with I/O control 
routine (F$IO) and standard magnetic tape routines. 

When F$IO is called, the format descriptor list and a flag 
indicating input are transmitted. 

When the F$IO routine needs a buffer of data, it calls this 
driver, which in turn calls the appropriate magnetic tape unit 
and conversion routines (for formatted READ). The number 
of words read is 60, equivalent to 120 characters. 

The appropriate magnetic tape units are physical magnetic tape 
units 1 through 5 corresponding respectively to logical tape 
units, numbers 5 through 9. 



Data Type 



Information is in ASCII if formatted, or in binary if 
unformatted. 



Other Routines 
Used 



F$IO, I$MA, I$MC, C$6T0£ 
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F$Rn 



Purpose 



To control the input drivers for variable input device numbers. 



DAP Calling 
Sequence 



LDA 
CALL 
DAC 
(Return) 



d 
F$Rn 



Location of device number 

Location of Format Descriptor List = 00000 
if format is binary 



FORTRAN 
References 



READ(x, f) list 
READ(x, f) 
READ(x) list 
READ(x) 



f = FORTRAN statement number; 

x = variable device number 1 through 9 

unformatted read where x = 1 , Z, 3, 5, 6, 7, 8, 9 



Method 



The value of d is checked for correct limits and is then used 
to determine the entry position of a Jump Table. The Jump 
Table transfers to the proper F$R subroutine. (Note that 
the entire F$R subroutine must be called into memory along 
with this subroutine, because there is no way of knowing 
in advance which drivers are required. ) 

If d does not equal a number from 1 to 9, the computer 

halts with a 1 in the A-register. The A-register may be changed 

manually to another device number; otherwise, the typewriter 

will be selected as the input device when START is pressed 

to continue processing. 

Other errors, such as parity, end of tape, etc. , cause the 
actions described in the appropriate F$R subroutine. 



Data Type 



Information is in ASCII if formatted, or in binary if 
unformatted. 



Other Routines 
Used 



F$R1, F$R2, F$R3, F$R5-9 
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Purpose To aid in debugging programs by printing the following: 

• Values of variables or array elements as they are 
being stored 

• Values of the DO parameters as they vary 

' • Locations of statement numbers as they are encountered 

• Content of IF statements as they are evaluated. 



DAP Calling CALL F$TR 

Sequence OCT The first three bits of the first argument 

OCT are: 

OCT 000 = = Statement number 

(Return) 001 = 1 = Integer 

010 = 2 = Real 

011 = 3 = Logical 

100 = 4 (Not Used) 

101 - S = Com til ex 

110 = 6 = Double precision 
The mode of an IF statement depends 
on the mode of the expression being 
evaluated. The remaining 13 bits of 
ARG1 and all the ARG2 and ARG3 are 
the name of the variable, the statement 
number, or blank. 



FORTRAN TRACE xl,x2, . . . xn Each x represents a variable or array 

rrarne ; 
TRACE n n represents a statement number 



References 



Method The A-register is stored in location zero (index register). 

The location of the argument is extracted and is interchanged 
with the index register, restoring the A-register. Sense 
Switch 4 is interrogated and if set, return is made to the 
calling program. If Sense Switch 4 is reset, trace mode is 
entered. The A- and B -registers are saved, and the buffer 
pointer is reset. The mode jump address is then set up. 
The 1- to 6-character name is then put in the format state- 
ment buffer. The return address for exit is calculated and 
stored. The 1- to 6-character name is interrogated to determine 
whether a statement number or a variable is being processed. 
If a variable name is being processed, an equal sign is stored 
in the format buffer following the name of the variable. An 
indirect jump is then taken, contingent upon whether the mode 
is real, integer, logical, complex, or double -precision. The 
following formats are moved to the format buffer, depending 
on the mode, to provide maximum representation of the 
variable. REAL G19>? 

INTEGER 17 

LOGICAL L2 

COMPLEX E15.7.H, E15. 7 

DOUBLE PRECISION D19. 11 
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F$TR cont. 



The output device and format location are provided to the 
appropriate I/O driver and communicated to F$IO. Then the 
variable name and equal sign, and its appropriate value, are 
printed. When an IF statement evaluation is being written 
the four characters IF( ) are printed, followed by an equal 
sign and the value of the expression in the parentheses. 

If a statement number is being printed, the number is 
bracketed by opening and closing parentheses, and is printed. 
Upon completion of each printout, the buffer is closed and 
the A- and B -registers restored. Return is made to the 
calling program. 



Data Type See calling sequence. 

of Arguments 



Other Routines F$W1, F$AR, F$CB, AC1, AC2, AC3 

Used 
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F$W1 



Purpose 



To control the typewriter output routine. 



DAP Calling 
Sequence 



CALL 

DAC 

(Return) 



F$W1 



Location of the format descriptor list 



FORTRAN 
References 



WRITE(1, f)list 
WRITE(l.f) 



f = FORTRAN statement number 



Method 



This subroutine connects the calling program with the I/O 
control subroutine (F$IO). Included in this subroutine is the 
driving logic needed to produce output on the typewriter. 

When F$IO is called, the location of the format descriptor 
list (if any), the entry location of the driver subroutine, and 
a flag indicating output are transferred. 

After the F$IO subroutine has generated a buffer full of data 
(72 characters), return is made to the driver output entry 
of this subroutine. At that time, the first character of the 
buffer is analyzed for proper line feed control as follows: 

Blank — Type a single carriage return 

followed by characters (2-72) 

— Issue two carriage returns, then 

type characters (2-72) 

1 — Skip to the top of the next page, 

then type characters (2-72) 

+ — No line advance, type characters (2-72) 

Others — Same as blank, except characters 

(1 -72) are typed 

This subroutine prints 60 lines per page and skips six lines 
to the top of the next page. The operator should start three 
lines down the first page in order to get correct spacing 
between pages. 



Data Type 



Only ASCII information is processed. 



Other Routines 
Used 



F$IO, 0$AC, 0$AF 
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Purpose 



To control the paper tape output routine. 



DAP Calling 
Sequence 



FORTRAN 
References 



CALL 


F$Vi 


DAC 


n 


(Return) 




WRITE(2, 


f) list 


WRITE(2, 


£) 


WRITE(2) 


list 



Location of format descriptor list 
(00000 if output mode is binary) 



f = FORTRAN statement number 
Unformatted WRITE to paper tape punch 



Method 



This subroutine connects the calling program to the I/O control 
subroutine F$IO. Included in the subroutine is the driving logic 
required to produce output on the paper tape punch. When F$IO 
is called, this subroutine transmits the location of the format 
descriptor list (if any) the entry location of the driver sub- 
routine, and a flag indicating input. 

After the F$IO subroutine has generated a full buffer of data 
60 words at two characters per word, or 40 words at three 
binary characters per word, return is made to the driver 
output entry of the subroutine. At that time, the buffer is 
punched on tape. 



Data Type 



Information is in ASCII if formatted or in binary if unformatted. 



Other Routines 
Used 



F$IO, 0$PF, 0$PP, 0$PB 
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Purpose 



To control the card punch routines. 



DAP Calling 
Sequence 



FORTRAN 
References 



CALL 

DAC 

(Return) 


F$W3 
n 


WRITE(3, 
WRITE(3, 
WRITE(3) 


f) list 

f) 
list 



Location of format descriptor list 
(00000 if output mode is binary) 



f = format statement number 
Unformatted WRITE to the card punch 



Method 



This subroutine connects the calling program (FORTRAN 

Object Program) to the I/O control subroutine, F$IO, and to 

the card punch subroutines. When F$IO is called, F$W3 

transmits the location of the address of the format descriptor 

list (if any), including a flag indicating output mode (DAC* for input), 

and a location for rpfntrance to F$W3. 

After the F$IO subroutine has generated a full buffer of data 
(40 words at two BCD characters per word or 60 binary words), 
return is made to F$W3. The appropriate card punch subroutine 
is called, and a card is punched. 



Data Type 



Information is in Hollerith if formatted, or in column binary 
if unformatted. 



Other Routines 
Used 



F$IO, 0$CH, 0$CB 



Note 



The Hollerith information can be either 026 or 029 character 
set, depending on the version of 0$CH selected. 
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F$W4 



Purpose 



To control the line printer output routine. 



DAP Calling 
Sequence 



CALL 

DAC 

(Return) 



F$W4 

n Location of format descriptor list 

If n is zero, which normally indicates binary- 
output, the computer halts. Push START to 
print data in BCD format. 



FORTRAN 
References 



WRITE (4, f) list 
WRITE (4, f) 



format statement number 



Method 



This subroutine connects the calling program to the I/O control 
subroutine F$IO. Included in the subroutine is the driving logic 
required to produce output on the line printer. When F$IO 
is called, this routine transmits the location of the format 
descriptor list (if any), the entry location of the driver sub- 
routine, and a flag indicating output mode. 

After the F$IO subroutine has generated a full buffer of data 
(120 characters), return is made to the driver output entry 
of this subroutine. At that time, the first character of the 
buffer is analyzed for proper line feed control as follows: 

Blank — Advance one line and print characters 

2 through 12 0. 

— Advance two lines and print characters 

2 through 120. 

1 — Advance to top of next page and print 

characters 1 through 120. 

+ — Print characters 1 through 120 without 

advancing line position. 

Others — Advance one line and print characters 

1 through 12 0. 



Data Type 



Information is in ASCII format. 



Other Routines 

Used 



F$IO, 0$LP, 0$LO 
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Pu r po s e 



To control writing on magnetic tape. 



DAP Calling 
Sequence 



CALL F$Wx x = 5, 6, 7, 8, or 9 

DAC n Address of format, if any 



FORTRAN 
References 



WRITE(x, f) list 
WRITE (x, f) 
WRITE(x) list 



f = FORTRAN statement number 

x = 5,6, 7, 8, or 9 

Unformatted READ, x = 5, 6, 7, 8, or 9 



Method 



This program connects the calling program with the I/O control 
routine (F$IO) and the standard magnetic tape routines. When 
F$IO is called, the format descriptor list and a flag indicating 
output mode are transferred. 

When the F$IO routine has a buffer of data to write, it calls 
the driver, which in turn calls the appropriate magnetic tape 



i ^ X «J J.W 1 



^ ^ u± - J HTDTTTTM 



60 words are written, equivalent to 120 characters in either 
mode (formatted or unformatted). 

The appropriate magnetic tape units (physical) are numbers 
1 through 5 corresponding respectively to the logical tape 
units numbers 5 through 9 given for x above. 



Data Type 



Information on the tape is in ASCII if formatted, or in 
binary if unformatted. 



Other Routines 
Used 



F$IO, 0$MC, 0$MA, C$8T06 
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F$Wn 



Purpose 



To control the output drivers for variable output device numbe 



DAP Calling 
Sequence 



FORTRAN 
References 



LDA 
CALL 
DAC 
(Return) 


d 

F$Wn 

n 


WRITE(x ; 
WRITE (x ; 
WRITE(x; 


, f) list 
,f) 
) list 



Location of device number 

Location of format descriptor list 
(00000 if format is binary) 



f = FORTRAN statement number 
x = variable device number 1 through 9 
Unformatted WRITE, where x = 1,2,3,4,5, 
6,7,8, or 9 



Method 



The value of d is checked for correct limits and then used to 
determine the entry position of a Jump Table. The Jump 
Table is then transferred to the proper F$W subroutine. 
(Note that all F$W subroutines must be called into memory 
along with this subroutine, because there is no way of knowing 
in advance which drivers are required. 

If d does not equal a number from 1 through 9, the computer 
halts with a 1 in the A-register. The A-register may be 
changed manually to another device number. Otherwise, 
the typewriter will be selected as the output device when 
START is pressed to continue processing. 

Other errors, such as parity, end of tape, etc. , cause the 
actions described in the appropriate F$W subroutine. 



Data Type 



Information is in ASCII if formatted, or in binary if unformatted. 



Other Routines 

Used 



F$W1, F$W2, F$W3, F$W4, F$W5-9 



6-2 5 



AM74 



APPENDIX A 
TAPE CONTENTS 

MAGNETIC TAPE 70182805000 - LTCSI S 
(LIBRARY SOURCES CODED IN FORTRAN) 

This tape consists of the individual sources of the following programs in the order listed. 
This tape is one of two distributed and contains that portion of the FORTRAN Library that is 
FORTRAN-coded. 



FILE NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 



NAME 

STMEAN 

STGEOM 

STCORR 

STSPER 

STCRMT 

STMEDT 

ST CHI 2 

STBNPB 

STL NRG 

STCHIS 

STPLRG 

STANVI 

STANV2 

STANVR 

STANVL 

STANVG 

STANVB 

STANVY 

DEFOAD 

DEFOMA 

DEFOHA 

DEFORK 

DESOAD 

DESOMA 

DESOHA 

DESORK 

PLYMUL 

PLYDIV 

PLYINT 

PLYIRT 

PLYEVL 

PLYDEF 

NACPLY 

NAAITK 

NALAGR 

NABAIR 

MATTRS 

MATMUL 

MARITH 



DOC. NO. 

70181386000 
70181387000 
70181388000 
70181389000 
70181390000 
70181391000 
70181392000 
70181394000 
70181395000 
70181393000 
70181396000 
70181397000 
70181398000 
70181399000 
70181400000 
70181401000 
70181402000 
70181403000 
70181405000 
70181404000 
70181406000 
70181407000 
70181408000 
70181409000 
70181410000 
70181411000 
70181416000 
70181417000 
70181420000 
70181421000 
70181422000 
70181423000 
70181429000 
70181418000 
70181419000 
70181424000 
70181412000 
70181413000 
70181414000 
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FILE NUMBER 



NAME 



DOC. NO. 



40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 



MATEIG 

MATINV 

NAREGU 

NAMULL 

GSEID 

SORT 

SORT2 

CVPOLR 

E$62 

E$61 

E$26 

E$66 

DSQRT 

DCOS 

DSIN 

DEXP 

DLOGIO 

DLOG 

DLOG2 

DATAN2 

DATAN 

DMOD 

DSIGN 

DABS 

A$62 

S$62 

M$62 

D$62 

C$16 

DBLE 

CSQRT 

CCOS 

CSIN 

CLOG 

CEXP 

CABS 

E$51 

A$52 

S$52 

M$52 

D$52 

A$55 

S$55 

M$55 

D$55 

CONJG 

C$25 

CMPLX 

N$55 

END 



70181415000 

70181427000 

70181425000 

70181426000 

70181428000 

70181430000 

70181431000 

70181432000 

70180053000 

70180052000 

70182582000 

70180054000 

70182580000 

70180055000 

70182583000 

70182581000 

70180051000 

70182579000 

70182914000 

70180056000 

70182584000 

/Ul«2b»»UUU 

70182589000 

70182587000 

70180037000 

70180038000 

70180039000 

70180040000 

70180059000 

70180058000 

70182592000 

70180066000 

70182595000 

70182591000 

70182593000 

70182596000 

70182594000 

70180041000 

70180042000 

70180043000 

70180044000 

70182544000 

70180093000 

70182545000 

70180034000 

70182598000 

70180068000 

70182597000 

70180069000 
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MAGNETIC TAPE 70182806000 - LTCS2S 
(LIBRARY SOURCES CODED IN DAP) 



FILE NUMBER 



NAME 



1 


DMAX1 


2 


DMIN1 


3 


DINT 


4 


Z$80 


5 


A$81 


6 


C$61 


7 


A$66 


8 


A$66XRA 


9 


H$66 


10 


C$26 


11 


H$55 


12 


MAXO 


13 


MAX1 


14 


MINO 


15 


MINI 


16 


TANH 


17 


SQRT 


18 


SQRTX 


19 


SIN, COS 


20 


ATAN 


21 


E$21 


22 


E$22 


23 


ALOG 


24 


ALOGX 


25 


EXP 


26 


E$ll 


27 


E$11X 


28 


ABS 


29 


C$62 


30 


AMOD 


31 


L$66 


32 


AINT 


33 


N$66 


34 


DIM 


35 


SIGN 


36 


AIMAG 


37 


L$55 


38 


IFIX 


39 


FLOAT 


40 


C$12 


41 


C$21 


42 


LOC 


43 


C$81 


44 


1ST ORE 


45 


N$33 


46 


IFETCH 


47 


IABS 


48 


F$OE 


49 


MOD 


50 


F$TR-RA 


51 


SUB$ 



DOC. NO. 

70182585000 
70182586000 
70180850000 
70180851000 
70180852000 
70182554000 
70180853000 
70180979000 
70180855000 
70180857000 
70180860000 
70182548000 
70182549000 
70180649000 
70182551000 
70182565000 
70182560000 
70180681000 
70182563000 
70182564000 
70182562000 
70180045000 
70182559000 
70180682000 
70182561000 
70182547000 
70180684000 
70182570000 
70180884000 
70182572000 
70180854000 
70182571000 
70180856000 
70182573000 
70182574000 
70180858000 
70180859000 
70182553000 
70180062000 
70182575000 
70182558000 
70181962000 
70180882000 
70181982000 
70180090000 
70181983000 
70182552000 
70181984000 
70182555000 
70180827000 
70185150000 
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FILE NUMBER 

5 2 
5 3 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 

74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 



NAME 


DOC. NO. 


FbGA 


7 I S 5 1 5 1000 


FSGC 


70J85152000 


IDIM 


70182556000 


A$22 


70182536000 


M$22 


70182537000 


A$22X11 


70181805000 


M$22X11 


70181806000 


D$22X11 


70181804000 


ISIGN 


70182557000 


L$22 


70182534000 


H$22 


70182535000 


N$22 


70180097000 


SLITE 


70182599000 


M$ll 


70180035000 


D$ll 


70182546000 


M$11X 


70180685000 


D$11X 


70180686000 


OVERFL 


70180894000 


F$AT 


70180071000 


L$33 


70180065000 


F$WN 


70180089000 


VtBN 


70180088000 


F$R1 


70182610000 


F$W1 


70182611000 


F$R2 


70182612000 


F$W2 


70182613000 


F$R3 


70182614000 


F$W3 


70181667000 


F$W4 


70182616000 


F$R5-9 


70180306000 


F$F5-9 


70180310000 


F$W5-9 


70180307000 


F$10 


70182618000 


016CHAIN 


70180659000 


ARG$ 


70180072000 


F$D5-9 


70180308000 


F$B5-9 


70180309000 


F$ER-RA 


70181068000 


AC1 


70180717000 


END 




E 70182803541 - 


- LTCM1S 



(LIBRARY OBJECTS - SOFTWARE VERSION) 

This magnetic tape consists of the concatenation of the individual objects of the listed 
programs. They have been translated by the FORTRAN Translator Mod 1 if FORTRAN coded 
and/or assembled by the DAP- 16 Mod 2 Assembler. 

DOC. NO. 



FILE NUMBER 

1 
2 
3 

4 
5 
6 
7 



NAME 

ST MEAN 

STGEOM 

STCORR 

STSPER 

STCRMT 

STMEDT 

STCH12 



70181386000 
70181387000 
70181388000 
70181389000 
70181390000 
70181391000 
70181392000 
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FILE NUMBER 



NAME 



DOC. NO. 



9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 



STBNPB 

STLNRG 

ST CHIS 

STPLRG 

STANV1 

STANV2 

STANVR 

STANVL 

STANVG 

STANVB 

STANVY 

DEFOAD 

DEFOMA 

DEFOHA 

DEFORK 

DESOAD 

DESOMA 

DESOHA 

DESORK 

PLYMUL 

PLYDIV 

PLYINT 

PLYIRT 

OLYEVL 

PLYDIF 

NACPLY 

NAAITK 

NALAGR 

NABAIR 

MATTRS 

MATMUL 

MARITH 

MATEIG 

MATINV 

NAREGU 

NAMULL 

GSEID 

SORT 

SORT 2 

CVPOLR 

E$62 

E$61 

E$26 

E$66 

DSQRT 

DCOS 

DSIN 

DEXP 

DLOG10 

DLOG 

DLOG2 

DATAN2 

DATAN 

DMOD 

DSIGN 

DABS 

A$62 



70181394000 
70181395000 
70181393000 
70181396000 
70181397000 
70181398000 
70181399000 
70181400000 
70181401000 
70181402000 
70181403000 
70181405000 
70181404000 
70181406000 
70181407000 
70181408000 
70181409000 
70181410000 
70181411000 
70181416000 
70181417000 
70181420000 
70181421000 
70181422000 
70181423000 
70181429000 
70181418000 
70181419000 
70181424000 
70181412000 
70181431000 
70181414000 
70181415000 
70181427000 
70181425000 
70181426000 
70181428000 
70181430000 
70181431000 
70181432000 
70180053000 
70180052000 
70182582000 
70180054000 
70182580000 
70180055000 
70182583000 
70182581000 
70180051000 
70182579000 
70182914000 
70180056000 
70182584000 
70182588000 
70182589000 
70182587000 
70180037000 
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TP MTTHTDT?D 



65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 



SS62 

MS62 

D$62 

C$16 

DBLE 

CSQRT 

CCOS 

CSIN 

CLOG 

CEXP 

CABS 

E$51 

A$52 

S$52 

M$52 

D$52 

A$55 

S$55 

M$55 

D$55 

CONJG 

C$25 

CMPLX 

N$55 

DMAX1 

DMIN1 

DINT 

Z$80 

A$81 

C$61 

A$66 

H$66 

C$26 

H$55 

MAXO 

MAX1 

MINO 

MINI 

TANH 

SQRT 

SIN, COS 

ATAN 

E$21 

E$22 

ALOG 

EXP 

E$ll 

ABS 

C$62 

AMOD 

L$66 

AINT 

N$66 

DIM 

SIGN 

AIMAG 



70180036000 

70180039000 

70180040000 

70180059000 

70180058000 

70182592000 

70180066000 

70182595000 

70182591000 

70182593000 

70182596000 

70182594000 

70180041000 

70180042000 

70180043000 

70180044000 

70182544000 

70180093000 

70182545000 

70180034000 

70182598000 

70180068000 

70182597000 

70180069000 

70182585000 

70182586000 

70180850000 

70180851000 

70180852000 

70182554000 

70180853000 

70180855000 

70180857000 

70180860000 

70182548000 

70182549000 

70180649000 

70182551000 

70182565000 

70182560000 

70182563000 

70182564000 

70182562000 

70180045000 

70182559000 

70182561000 

70182547000 

70182570000 

70180884000 

70182572000 

70180854000 

70182571000 

70180856000 

70182573000 

70182574000 

70180858000 
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FILE NUMBER 


NAME 


DOC. NO. 


121 


L$55 


70180859000 


122 


IF IX 


70182553000 


123 


FLOAT 


70180062000 


124 


C$12 


70182575000 


125 


C$21 


70182558000 


126 


LOC 


70181962000 


127 


C$81 


70180882000 


128 


ISTORE 


70181982000 


129 


N$33 


70180090000 


130 


IFETCH 


70181983000 


131 


IABS 


70182552000 


132 


F$DE 


70181984000 


133 


MOD 


70182555000 


134 


F$TR-RA 


70180827000 


135 


SUB$ 


70185150000 


136 


F$GA 


70185151000 


137 


F$GC 


70185152000 


138 


IDIM 


70182556000 


139 


A$22 


70182536000 


140 


M$22 


70182537000 


141 


ISIGN 


70182557000 


142 


L$22 


70182534000 


143 


H$22 


70182535000 


144 


N$22 


70180097000 


145 


SLITE 


70182599000 


146 


M$ll 


70180035000 


147 


D$ll 


70182546000 


148 


OVERFL 


70180894000 


149 


F$AT 


70180071000 


150 


L$33 


70180065000 


151 


F$WN 


70180089000 


152 


F$RN 


70180088000 


153 


F$R1 


70182610000 


154 


F$W1 


70182611000 


155 


F$R2 


70182612000 


156 


F$W2 


70182613000 


157 


F$R3 


70182614000 


158 


F$W3 


70181667000 


159 


F$W4 


70182616000 


160 


F$R5-9 


70180306000 


161 


F$F5-9 


70180310000 


162 


F$W5-9 


70180307000 


163 


F$10 


70182618000 


164 


016CHAIN 


70180659000 


165 


ARG$ 


70180072000 


166 


F$D5-9 


70180308000 


167 


F$B5-9 


70180309000 


168 


F$ER-RA 


70181068000 


169 


AC1 


70180717000 


170 


SQRX1 


70188775000 


171 


COSX1 


70188781000 


172 


SINX1 


70188777000 


173 


ANT XI 


70188779000 


174 


LGEX1 


70188814000 


175 


LG2X1 


70188784000 


176 


EXEX1 


70188786000 


177 


EX2X1 


70188782000 



A-7 



AM 74 



FILE NUMBER 



178 


179 


180 


181 


182 


183 


184 


185 


186 


187 


188 


189 


190 


191 


192 


193 


Files 


Files 


£ lies 



NAME 

DSQRX1 

DCOSX1 

DSINX1 

DAT NX 1 

DLGEX1 

DLG2X1 

DEXEX1 

DEX2X1 

DMPY 

MPY 

DIV 

DADD 

DSUB 

ROND 

RODD 

TWOS 

END 

1-47 
48-169 

1 (U-LV3 



DOC. NO. 

70188788000 
70188792000 
70188790000 
70188793000 
70188801000 
70188795000 
70188799000 
70188797000 
70188808000 
70188811000 
70188810000 
70188812000 
70188813000 
70188805000 
70188804000 
70188803000 



Statistical Library 
FORTRAN Library 
£ ixed J^omt iviath Libra. 



MAGNETIC TAPE 70182804541 - LTCM1H 
(LIBRARY OBJECTS - HARDWARE VERSION) 

This magnetic tape consists o£ the concatenation of the individual objects of the listed 
programs. They have been translated by the FORTRAN Translator Mod 1 if FORTRAN cod< 
and/or assembled by the DAP-16 Mod 2 Assembler. 



FILE NUMBER 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 



NAME 

STMEAN 

STGEON 

STCORR 

STSPER 

STCRMT 

STMEDT 

STCHI2 

STBNPB 

STLNRG 

STCHIS 

STPLRG 

STANV1 

STANV2 

STANVR 

STANVL 

STANVG 

STANVB 

STANVY 

DEFOAD 

DEFOMA 

DEFOHA 

DEFORK 

DESOAD 



DOC. NO. 

70181386000 
70181387000 
70181388000 
70181389000 
70181390000 
70181391000 
70181392000 
70181394000 
70181395000 
70181393000 
70181396000 
70181397000 
70181398000 
70181399000 
70181400000 
70181401000 
70181402000 
70181403000 
70181405000 
70181404000 
70181406000 
70181407000 
70181408000 
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FILE NUMBER 



NAME 



DOC. NO. 



24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 



DESONA 

DESOHA 

DESORK 

PYLMUL 

PLYDIV 

PLYINT 

PLYIRT 

PLYEVL 

PLYDIF 

NACPLY 

NAAITK 

NALAGR 

NAB AIR 

NATTRS 

NATMUL 

MARITH 

NATEIG 

NATENV 

NAREGU 

NAMULL 

GSEID 

SORT 

SORT 2 

CVPOLR 

E$62 

E$61 

E$26 

E$66 

DSQRT 

DCOS 

DSIN 

DEXP 

DLOG10 

DLOG 

DLOG2 

DATAN2 

DATAN 

DMOD 

DSIGN 

DABS 

A$62 

S$62 

M$62 

D$62 

C$16 

DBLE 

CSQRT 

CCOS 

CSIN 

CLOG 

CEXP 

CABS 

E$51 

A$52 

S$52 

M$52 

D$52 



70181409000 

70181410000 

70181411000 

70181416000 

70181417000 

70181420000 

70181421000 

70181422000 

70181423000 

70181429000 

70181418000 

70181419000 

70181424000 

70181412000 

70181413000 

70181414000 

70181415000 

70181427000 

70181425000 

70181426000 

70181428000 

70181430000 

70181431000 

70181432000 

70180053000 

70180052000 

70182582000 

70180054000 

70182580000 

70180055000 

70182583000 

70182581000 

70180051000 

70182579000 

70182914000 

70180056000 

70182584000 

70182588000 

70182589000 

70182587000 

70180037000 

70180038000 

70180039000 

70180040000 

70180059000 

70180058000 

70182592000 

70180066000 

70182595000 

70182591000 

70182593000 

70182596000 

70182594000 

70180041000 

70180042000 

70180043000 

70180044000 



A- 9 
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TPTT T7> TVTTTA A O VT D 

81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 



NAME 


DOC. NO. 


A$55 


70182544000 


S$55 


70180093000 


M$55 


70182545000 


D$55 


70180034000 


CONJG 


70182598000 


C$25 


70180068000 


CMPLX 


70182597000 


N$55 


70180069000 


DMAX1 


70182585000 


DMIN1 


70182586000 


DINT 


70180850000 


Z$80 


70180851000 


A$81 


70180852000 


C$61 


70182554000 


A$66XKA 


70180979000 


H$66 


70180855000 


C$26 


70180857000 


H$55 


70180860000 


MAXO 


70182548000 


MAX1 


70182549000 


MING 


70180649000 


MINI 


70182551000 


TANH 


70182565000 


SQRTX 


70180681000 


SIN, COS 


70182560000 


ATAN 


70182564000 


E$21 


70182562000 


E$22 


70180045000 


ALOGX 


70180682000 


EXP 


70182561000 


E$11X 


70180684000 


ABS 


70182570000 


C$62 


70180884000 


AMOD 


70182572000 


L$66 


70180854000 


A INT 


70182571000 


N$66 


70180856000 


DIM 


70182573000 


SIGN 


70182874000 


AIMAG 


70180858000 


L$55 


70180859000 


IF DC 


70182553000 


FLOAT 


70180062000 


C$12 


70182575000 


C$21 


70182558000 


LOC 


70181962000 


C$81 


70180882000 


ISTORE 


70181982000 


N$33 


70180090000 


IFETCH 


70181983000 


IABS 


70182552000 


F$OE 


70181984000 


MOD 


70182555000 


F$TR-RA 


70180827000 


SUB$ 


70185150000 


F$GA 


70185151000 


F$GC 


70185152000 
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ilJ_.ii; WUMBER 


NAME 


DOC. NO. 


138 


IDIM 


70182556000 


139 


A$22X11 


70181805000 


140 


M$22X11 


70181806000 


141 


D$22X11 


70181804000 


142 


ISIGN 


70182557000 


143 


L$22 


70182534000 


144 


H$22 


70182535000 


145 


N$22 


70180097000 


146 


SLITE 


70182599000 


147 


M$11X 


70180685000 


148 


D$11X 


70180686000 


149 


OVERFL 


70180894000 


150 


F$AT 


70180071000 


151 


L$33 


70180065000 


152 


F$WN 


70180089000 


153 


F$RN 


70180088000 


154 


F$R1 


70182610000 


155 


F$W1 


70182611000 


156 


F$R2 


70182612000 


157 


F$W2 


70182613000 


158 


F$R3 


70182614000 


159 


F$W3 


70181667000 


160 


F$W4 


70182616000 


161 


F$R5-9 


70180306000 


162 


F$F5-9 


70180310000 


163 


F$W5-9 


70180307000 


164 


F$I0 


70182618000 


165 


016CHAIN 


70180659000 


166 


ARG$ 


70180072000 


167 


F$D5-9 


70180308000 


168 


F$B5-9 


70180309000 


169 


F$ER-RA 


70181068000 


170 


AC1 


70180717000 


171 


SQRX2 


70188776000 


172 


COSX2 


70180761000 


173 


SINX2 


70188778000 


174 


ATNX2 


70188780000 


175 


LGEX2 


70188815000 


176 


LG2X2 


70188785000 


177 


EXEX2 


70188787000 


178 


EX2X2 


70188783000 


179 


DSQRX2 


70188789000 


180 


DCOSX2 


70180762000 


181 


DSINX2 


70188791000 


182 


DATNX2 


70188794000 


183 


DLGEX2 


70188802000 


184 


DLG2X2 


70188796000 


185 


DEXEX2 


70188800000 


186 


DEX2X2 


70188798000 


187 


DMPYH 


70188809000 


188 


DADD 


70188812000 


189 


DSUB 


70188813000 


190 


ROND 


70188805000 


191 


RODD 


70188804000 


192 


TWOS 
END 


70188803000 
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Files 
Files 
Files 



1-47 

48-170 

171-192 



Statistical Library 
FORTRAN Library 
Fixed Point Math Library 



PAPER TAPE 70181876000 - LTCF1 (Paper Tape 1 of 6) 
This paper tape consists of the concatenation of the individual objects of the listed 
programs. They have been translated by the FORTRAN Translator MOD 1 and assembled by 
the DAP- 16 Mod 2 Assembler. 



FILE NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 



NAME 

E$62 

E$61 

E$26 

E$66 

DSQRT 

DCOS 

DSIN 

DEXP 

DLOG10 

DLOG 

DLOG2 

DATAN2 

DATAN 

DMOD 

DSIGN 

DABS 

A$62 

S$62 

M$62 

D$62 

C$16 

DBLE 

END 



DOC. NO. 

70180053000 

70180052000 

70182582000 

70180054000 

70182580000 

70180055000 

70 182583000 

70182581000 

70180051000 

70182579000 

70182914000 

70180056000 

70182584000 

70182588000 

70182589000 

70182587000 

70180037000 

70180038000 

70180039000 

70180040000 

70180059000 

70180058000 



PAPER TAPE 70181877000 - LTCF2 (Paper Tape 2 of 6) 
This paper tape consists of the concatenation of the individual objects of the listed 
program. They have been translated by the FORTRAN Translator MOD 1 and assembled by 
the DAP-16 Mod 2 Assembler. 



FILE NUMBER 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



NAME 

CSQRT 

CCOS 

CSIN 

CLOG 

CEXP 

CABS 

E$51 

A$52 

S$52 

M$52 



DOC. NO. 

70182592000 
70180066000 
70182595000 
70182591000 
70182593000 
70182596000 
70182594000 
70180041000 
70180042000 
70180043000 
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FILE NUMBER 



11 


12 


13 


14 


15 


16 


17 


18 


19 



NAME 

D$52 

A$55 

S$55 

M$55 

D$55 

CONJG 

C$25 

CMPLX 

N$55 

END 



DOC. NO. 

70180044000 
70182544000 
70180093000 
70182545000 
70180034000 
70182598000 
70180068000 
70182597000 
70180069000 



PAPER TAPE 70181882000 - LTCF3S (Paper Tape 3 of 6) 
This is the Software Version of Paper Tape 3 



FILE NUMBER 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 



NAME 

DMAX1 

DMIN1 

DINT 

Z$80 

A$81 

C$61 

A$66 

H$66 

C$26 

H$55 

MAX0 

MAX1 

MIN0 

MINI 

TANH 

SQRT 

SIN, COS 

ATAN 

E$21 

E$22 

ALOG 

EXP 

E$ll 

END 



DOC. NO. 

70182585000 

70182586000 

70180850000 

70180851000 

70180852000 

70182554000 

70180853000 

70180855000 

70180857000 

70180860000 

70182548000 

70182549000 

70180649000 

70182551000 

70182565000 

70182560000 

70182563000 

70182564000 

70182562000 

70180045000 

70182559000 

70182561000 

70182547000 



PAPER TAPE 70181878000 - LTCF3H (Paper Tape 3 of 6) 
This is the Hardware Version of Paper Tape 3 



FILE NUMBER 



NAME 



1 


DMAX1 


2 


DMIN1 


3 


DINT 


4 


Z$80 


5 


A$81 


6 


C$61 


7 


A$66XRA 


8 


H$66 


9 


C$26 



DOC. NO. 

70182585000 
70182586000 
70180850000 
70180851000 
70180852000 
70182554000 
70180979000 
70180855000 
70180857000 
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FILE NUMBER 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 



NAME 

H$55 

MAXO 

MAX1 

MINO 

MINI 

TANH 

SQRTX 

SIN, COS 

ATAN 

E$21 

E$22 

ALOGX 

EXP 

E$11X 

END 



DOC. NO. 

70180860000 
70182548000 
70182549000 
70180649000 
70182551000 
70182565000 
70180681000 
70182563000 
70182564000 
70182562000 
70180045000 
70180682000 
70182561000 
70180684000 



PAPER TAPE 70181879000 - LTCF4 (Paper Tape 4 of 6) 



FILE NUMBER 

1 
2 
3 
4 
5 
6 
7 



9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



NAME 
ABS 

AMOD 

L$66 

AINT 

N$66 

DIM 

SIGN 

AIMAG 

L$55 

IFIX 

FLOAT 

C$12 

C$21 

LOC 

C$81 

ISTORE 

N$33 

IF ETCH 

LABS 

F$OE 

MOD 

F$TR-RA 

SUB$ 

F$GA 

F$GC 

END 



DOC. NO. 

70182570000 

70 180884000 

70182572000 

70180854000 

70182571000 

70180856000 

70182573000 

70182574000 

70180858000 

70180859000 

70182553000 

70180062000 

70182575000 

70182558000 

70181962000 

70180882000 

70181982000 

70180090000 

70181983000 

70182552000 

70181984000 

70182555000 

70180827000 

70185150000 

70185151000 

70185152000 



PAPER TAPE 70181883000 - LTCF5S (Paper Tape 5 of 6) 
This is the Software Version of Tape 5 

DOC. NO. 



FILE NUMBER 

1 
2 
3 

4 



NAME 

IDIM 
A$22 
M$22 
ISIGN 



70182556000 
70182536000 
70182537000 
70182557000 
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FILE NUMBER 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



NAME 

L$22 

H$22 

N$22 

SLITE 

M$ll 

D$ll 

OVERFL 

F$AT 

L$33 

F$WN 

F$RN 

F$R1 

F$W1 

F$R2 

F$W2 

F$R3 

F$W3 

F$W4 

F$R5-9 

F$F5-9 

F$W5-9 

F$IO 

END 



DOC. NO. 

70182534000 

70182535000 

70180097000 

70182599000 

70180035000 

70182546000 

70180894000 

70180071000 

70180065000 

70180089000 

70180088000 

70182610000 

70182611000 

70182612000 

70182613000 

70182614000 

70181667000 

70182616000 

70180306000 

70180310000 

70180307000 

70182618000 



PAPER TAPE 70181880000 - LTCF5H fPaper Tape 5 of 6) 
This is the Hardware Version of Tape 5 



FILE NUMBER 



NAME 



1 


IDIM 


2 


A$22X11 


3 


M$22X11 


4 


D$22X11 


5 


ISIGN 


6 


L$22 


7 


H$22 


8 


N$22 


9 


SLITE 


10 


M$11X 


11 


D$11X 


12 


OVERFL 


13 


F$AT 


14 


L$33 


15 


F$WN 


16 


F$RN 


17 


F$R1 


18 


F$W1 


19 


F$R2 


20 


F$W2 


21 


F$R3 


22 


F$W3 


23 


F$W4 


24 


F$R5-9 


25 


F$F5-9 


26 


F$W5-9 


27 


F$IO 




END 



DOC. NO. 

70182556000 

70181805000 

70181806000 

70181804000 

70182557000 

70182534000 

70182535000 

70180097000 

70182599000 

70180685000 

70180686000 

70180894000 

70180071000 

70180065000 

70180089000 

70180088000 

70182610000 

70182611000 

70182612000 

70182613000 

70182614000 

70181667000 

70182616000 

70180306000 

70180310000 

70180307000 

70182618000 



A-15 



AM 74 



PAPER TAPE 70181881000 - LTCF6 (Paper Tape 6 of 6) 



FILE NUMBER 


NAME 


DOC. NO. 


1 


016CHA1N 


70180659000 


2 


ARG$ 


70180072000 


3 


F$D5-9 


70180308000 


4 


F$B5-9 


70180309000 


5 


F$ER-RA 


70181068000 


6 


ACl 

END 


70180717000 
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APPENDIX B 
MATHEMATICAL ROUTINES 



Function 



Complex 



Absolute value 

Add 

Add real argument 

Conjugate 

Convert imaginary part to real 

Cosine 

Divide 

Divide by real argument 

Exponential, base e 

Load 

Load real part 

Logarithm, base e 

Multiply 

Multiply by real argument 

Negate 

Raise to integer power 

Sine 

Square root 

Store (hold) 

Subtract 

Subtract single -precision argument 



Routine 



CABS 

A$55 

A$52 

CONJG 

AIMAG 

CCOS 

D$55 

D$52 

CEXP 

L$55 

REAL 

CLOG 

M$55 

M$52 

N$55 

E$51 

CSIN 

CSQRT 

H$55 

S$55 

S$52 



Double -Precis ion 



Absolute value 

Add 

Add single-precision argument 

Add integer to exponent 

Arctangent, principal value 



DABS 

A$66 

A$62 

A$81 

DATAN 



B-l 
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Function 



Routine 



Double -Precision 



Arctangent, X/Y 

Clear (zero) exponent 

Convert exponent to integer 

Convert to integer 

Convert to single-precision 

Cos ine 

Divide 

Divide by real argument 

Exponential, base e 

Load 

Logarithm, base e 

Logarithm, base 2 

Logarithm, base 10 

Minimum value 
Multiply- 
Multiply by real argument 
Negate 

Raise to double-precision power 
Raise to integer power 
Raise to real power 
Remainder 
Sine 

Square root 
Store (hold) 
Subtract 

Subtract real argument 

Transfer sign of second argument to first 
Truncate fractional bits 

Truncate fractional bits and convert 
to integer 



DATAN2 

Z$80 

C$81 

C$61 

C$62 

DCOS 

D$66 

D$62 

DEXP 

L$66 

DLOG 

DLOG2 

DLOG10 

DMAX 1 

DMIN1 

M$66 

M$62 

N$66 

E$66 

E$61 

E$62 

DMOD 

DSIN 

DSQRT 

H$66 

S$66 

S$62 

DSIGN 

DINT 

ID INT 



Real 



Absolute value 

Add 

Arctangent, principal value 

Arctangent, X/Y 

Convert pair to complex 

B-2 



ABS 

A$22 
ATAN 
ATAN2 
CMPLX 



AM74 



Function 



Routine 



Real 










Convert to complex format 


C$25 




Convert to double-precision 


C$26 




Convert to integer 




C$21 




Divide 




D$22 




Exponential, base e 




EXP 




Hyperbolic tangent 




TANH 




Load 




L$22 




Logarithm, base e 




ALOG 




Logarithm, base 10 




ALOG10 




Maximum integer value 




MAX1 




Maximum value 




AMAX1 




Minimum integer value 




MINI 




Minimum value 




AMIN1 




Multiply- 




M$22 




Positive difference 




DIM 




Raise to double-precision power 


E$26 




Raise to integer power 




E$21 




Raise to real power 




E$22 




Remainder 




AMOD 




Sine, cosine 




SIN, COS 




Square root 




SQRT 




Store (hold) 




H$22 




Subtract 




S$22 




Transfer sign of second 


argument to first 


SIGN 




Truncate fractional bits 




A INT 




Truncate fractional bits 


and convert 






to integer 




IF IX, INT 




TWOs complement 




N$22 



Integer 



Absolute value 

Convert to double -precis ion 

Convert to real (FORTRAN-generated) 

Convert to real 

Divide 

Maximum value 



IABS 

C$16 

FLOAT 

C$12 

D$ll 

AMAXO 



B-3 



AM 74 



Function 



Integer 



Maximum integer value 

Minimum value 

Minimum integer value 

Multiply 

Positive difference 

Raise to integer power 

Remainder 

Transfer sign of second argument to first 



Routine 



MAXO 

AMINO 

MINO 

M$ll 

IDIM 

E$ll 

MOD 

ISIGN 



Logical 



Complement 

OR with A-register 



N$33 
L$33 



B-4 



AM74 



APPENDIX C 
SUBROUTINE FUNCTIONS 



INTRINSIC AND EXTERNAL FUNCTIONS 



Mathematical and Trigonometric Functions 



Name 

SIN 

DSIN 

CSIN 

COS 

DCOS 

CCOS 

AT AN 
DATAN 
ATAN2 
DATAN2 

TANH 

SQRT 

DSQRT 

CSQRT 

EXP 

DEXP 

CEXP 

ALOG 
DLOG 
CLOG 

ALOG 10 

DLOG2 

DLOG10 

ABS 
IABS 
DABS 
CABS 

AMOD 

MOD 

DMOD 

AINT 

DINT 

IDINT 

IFIX 

INT 



Argument 


Result 


Data 


Type 


Data Type 


R 




R 


D 




D 


C 




C 


R 




R 


D 




D 


C 




C 


R 




R 


D 




D 


R 




R 


D 




D 


R 




R 


R 




R 


D 




D 


C 




C 


R 




R 


D 




D 


C 




C 


R 




R 


D 




D 


C 




C 


R 




R 


D 




D 


D 




D 


R 




R 


I 




I 


D 




D 


C 




R 


R 




R 


I 




I 


D 




D 


R 




R 


D 




I 


D 




I 


R 




I 


R 




I 



Function 
Sine (radians) 

Cosine (radians) 

Arctangent (radians) 

Hyperbolic tangent (radians) 
Square root 

Exponential 

Natural logarithm 

Common logarithm 

Absolute value 

Remainder 

Truncate fractional bits 



C-l 



AM74 



Name 

AMAXO 

AMAX1 

MAXO 

MAX1 

DMAX1 

AMINO 
AMIN1 
DMIN1 
MIN10 
MINI 

FLOAT 

AIMAG 

DBLE 

CMPLX 

REAL 

SNGL 

SIGN 
DSIGN 

1SIGN 

DIM 
IDIM 

CONJG 



Argument 


R. 


esult 


Data Type 


D; 


ata Type 


I 




R 


R 




R 


I 




I 


R 




I 


D 




D 


I 




R 


R 




R 


D 




D 


I 




I 


R 




I 


I 




R 


C 




R 


R 




D 


C 




R 


C 




R 


R 




D 


R 




R 


D 




D 


I 




I 


R 




R 


I 




I 



Function 



Choose largest argument 



Choose smallest argument 



Change data type or argument 



Value of first argument, sign of 
second 

Positive difference 
Complex conjugate 



Special Subroutines for FORTRAN Use 



Name 

IFETCH(I) 

ISTORE(I, J) 

LOC 

OVERFL 

SLITE 

SLITET 

SSWTCH 



Argument 
Data Type 



Result 
Data Type 



Function 

Get contents of location I 
Store value of J in location I 
Find address of argument 
Check for error condition 
Set and reset sense lights or 
switches 



COMPILER SUPPORT SUBROUTINES 



Conversion Routines 



Name 



C$12 
C$16 
C$21 
C$25 
C$26 
C$61 
C$62 
C$81 



Argument 


Result 


Data 


Type 


Data Type 


I 




R 


I 




D 


R 




I 


R 




C 


R 




D 


D 




I 


D 




R 


D 




D 



Function 

Convert integer to real 
Convert integer to double -precision 
Convert real to integer 
Convert real to complex 
Convert real to double -precision 
Convert double -precision to integer 
Convert double -precision to real 
Convert exponent of double - 
precision number to integer 



C-2 



AM74 



Arithmetic Routines 



Name 


Function 


A$22 


R = R+R 


A$52 


C = C+R 


A$55 


C = C+C 


A$62 


D = D+R 


A$66 


D = D+D 


A$81 


D = D*(2**I) 


D$ll 


1 = 1/1 


D$22 


R = R/R 


D$52 


C = C/R 


D$55 


C = C/C 


D$62 


D = D/R 


D$66 


D = D/D 


E$ll 


I = 1**1 


E$21 


R = R**I 


E$22 


R = R**R 


E$26 


D = R**D 


E$51 


C = C**I 


E$61 


D = D**I 



Name 


Function 


E$62 


D = D**R 


E$66 


D = D**D 


M$ll 


I = 1*1 


M$22 


R = R*R 


M$52 


C = C*R 


M$55 


C = C*C 


M$62 


D = D*R 


M$66 


D = D*D 


N$22 


R = -R 


N$33 


L = -L 


N$55 


C = -C 


N$66 


D = -D 


S$22 


R = R-R 


S$52 


C = C-R 


S$55 


C = C-C 


S$62 


D = D-R 


S$66 


D = D-D 


Z$80 


Replace binary exponent 




with zero 



Miscellaneous Routines 



Name 
AC1 

ARG$ 

H$22 
H$55 
H$66 

L$22 
L$55 
L$66 

L$33 

SUB$ 



Function 

Pseudo accumulators 

Convert indirect address to direct address 

Store real number in memory- 
Store complex number in memory 
Store double -precision number in memory 

Load real number into A- and B -registers 

Load complex number into complex accumulator 

Load double -precision number into double-precision accumulator 

INCLUSIVE OR with A-register 

Calculate address of array element 



C-3 



AM74 



APPENDIX D 
LIBRARY INDEX 











Approx. 




Primary 


Entry 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words ) 
150 


Numbe r 


A$22 


A$22 


ARG$ 




5S 




S$22 


N$22 
F$ER 








A$22X 


A$22 
S$22 


N$22 
F$ER 




140 


5H 


A$52 


A$52 


F$AT 

H$55 

L$22 

A$22 

H$22 

L$55 




20 


2 


A$55 


A$55 


F$AT 

H$55 

SUB$ 

L$22 

A$22 

H$22 

L$55 


4 
2 
2 
2 
1 


60 


2 


A$62 


A$62 


F$AT 
H$66 
DBLE 
A$66 




20 


1 


A$66 


A$66 
S$66 
M$66 
D$66 


N$66 

F$ER 

H$66 

L$66 

ARG$ 

ACl 

AC2 

AC3 


11 


580 


3S 


A$66X 


A$66 

A$66X 

S$66 

S$66X 

M$66 

M$66X 

D$66 

D$66X 


N$66 
F$ER 
H$66 
L$66 
ARG$ 
ACl 
AC 2 
AC 3 


11 


530 


3H 



D-l 



AM74 











Approx. 




Primary 


Entry 


Subroutines 


Number 


of Storage 


Tape 


Name 


Points 


Called 


References (Words ) 


Number 


A$81 


A$81 


N$22 
F$ER 
AC1 
AC 2 


2 


70 


3 


ABS 


ABS 


L$22 

N$22 




10 


4 


AC1 


AC1 
AC2 
AC 3 
AC4 
AC5 






5 


6 


AIMAG 


AIMAG 


L$55 
L$22 
AC 3 




10 


4 


AINT 


AINT 


L$22 
N$22 
A$22 
S$22 




30 


4 


A LOG 


ALOG10 
A LOG 


ARG$ 

C$12 

H$22 

L$22 

A$22 

S$22 

D$22 

M$22 

F$ER 


5 
3 
6 
2 
1 
7 
1 


120 


3S 


ALOGX 


ALOG10 


ARG$ 


1 


180 


3H 




ALOG 


C$12 


1 








ALOGX 


A$22 
M$22 
S$22 

F$ER 


4 
4 
1 
1 






ALOG10 


See ALOG 


or ALOGX 








AMAXO 


See MAXO 










AMAX1 


See MAX1 










AMINO 


See MINO 










AMIN1 


See MINI 










A MOD 


AMOD 


L$22 
D$22 
AINT 
M$22 
N$22 
A$22 




30 


4 



ARG$ 



ARG$ 



20 



D-2 



AM74 











Approx. 




Primary 


Entry 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words 1Q ) 


Number 


ATAN 


ATAN2 


ARG$ 


3 


340 


3 




ATAN 


D$22 
N$22 
M$22 
A$22 
S$22 


6 
7 
5 
11 
2 






ATAN2 


See ATAN 










C$12 


C$12 


A$22 
N$22 




30 


4 


C$16 


C$16 


C$12 
C$26 




5 


1 


C$21 


C$21 


N$22 
A$22 
F$ER 




30 


4 


C$25 


C$25 


H$22 
CMPLX 




20 


2 


C$26 


C$26 


ACl 
AC2 
AC 3 




10 


3 


C$61 


C$61 


C$62 
C$21 




4 


3 


C$62 


C$62 
SNGL 


L$22 

N$66 

N$22 

L$66 

ACl 

AC2 




20 


4 


C$81 


C$81 


ACl 




10 


4 


CABS 


CABS 


F$AT 

SUB$ 

L$22 

M$22 

H$22 

A$22 

SQRT 


2 
2 
2 
2 


40 


2 


CCOS 


CCOS 


F$AT 

L$55 

A$55 

H$55 

CSIN 




40 


2 


CEXP 


CEXP 


F$AT 

SUB$ 

EXP 

H$22 

COS 

M$22 

SIN 

L$55 


2 


60 


2 



D-3 



AM 74 











Approx. 




Primary 


Entry 


Subroutines 


Number of 


Storage 


1 ape 


Name 


Points 


Called 


References 


(Words 1Q ) 
90 


Number 


CLOG 


CLOG 


F$AT 


1 


2 






SUB$ 


6 










L$22 


3 










M$22 


3 










H$22 


5 










A$22 


1 










ALOG 


1 










ATAN2 


1 










L$55 


1 






CMPLX 


CMPLX 


F$AT 
SUB$ 
L$22 
H$22 

L$55 


1 
2 
2 
2 

1 


40 


2 


CONJG 


CONJG 


F$AT 
SUB$ 
L$22 

H$22 
N$22 


1 

4 
2 
2 

1 


40 


2 


COS 


See SIN 


L$55 








CSIN 


CSIN 


F$AT 

SUB$ 

EXP 

H$22 

L$22 

D$22 

A$22 

SIN 

M$22 

S$22 

COS 

L$55 


1 
5 
1 
6 
3 
1 
1 
1 
4 
1 
1 
1 


90 


2 


CSQRT 


CSQRT 


F$AT 

SUB$ 

CABS 

H$22 

ABS 

A$22 

M$22 

SQRT 

L$22 

D$22 

L$55 


1 
7 
1 
8 
1 
1 
2 
1 
6 
1 
1 


90 


2 


D$ll 


D$ll 


ARG$ 
F$ER 


1 
1 


80 


5S 


D$11X 


D$ll 


ARG$ 


1 


40 


5H 




D$11X 


F$ER 


1 






D$22 


See M$22 










D$22X 


D$22 


N$22 
F$ER 


3 
2 


110 


5H 



D-4 



AM74 











Approx. 




Primary- 


Entry- 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words ) 
50 


Number 


D$52 


D$52 


F$AT 


1 


2 






H$55 


1 










SUB$ 


2 










L$22 


2 










D$22 


2 










H$22 


2 










L$55 


1 






D$55 


D$55 


F$AT 

H$55 

SUB$ 

L$22 

M$22 

H$22 

A$22 

D$22 

S$22 

N$22 

L$55 


1 
1 
12 
8 
6 
8 
2 
2 
1 
1 
1 


140 


2 


D$62 


D$62 


F$AT 

H$66 

DBLE 

L$66 

D$66 


1 
2 
1 
1 

1 


20 


1 


D$66 


See A$66 










D$66X 


See A$66X 










DABS 


DABS 


F$AT 

L$66 

N$66 




10 


1 


DATAN 


DATAN 


F$AT 

DABS 

H$66 

C$81 

L$66 

A$66 

N$66 

D$66 

M$66 


9 
1 
13 
10 
3 
2 
9 


180 


1 


DATAN2 


DATAN2 


F$AT 

L$66 

H$66 

F$ER 

D$66 

DATAN 

S$66 

A$66 


1 
9 
3 


70 


1 


DBLE 


DBLE 


F$AT 

L$22 

C$26 




20 


1 



D-5 



AM 74 











Approx. 




Primary 


Entry 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words 1Q ) 


Number 


DCOS 


DCOS 


F$AT 
L$66 
A$66 
H$66 
DSIN 


1 
1 
1 
1 
1 


20 


1 


DEXP 


DEXP 


F$AT 

L$66 

M$66 

H$66 

C$61 

C$16 

N$66 

A$66 

S$66 

D$66 

A$81 


1 
12 

8 
11 
1 
1 
1 
8 
3 
1 
1 


160 


1 


DIM 


DIM 


L$22 
S$22 


1 
1 


20 


4 


DINT 


DINT 


L$66 
N$66 
A$66 
S$66 
ACl 


1 
2 
1 
1 
1 


20 


3 


DIV$ 


See M$22 










DLOG 


DLOG 


F$AT 

DLOG2 

M$66 


1 

1 
1 


10 


1 


DLOG2 


DLOG2 


F$AT 

L$66 

F$ER 

C$81 

C$16 

H$66 

Z$80 

A$66 

S$66 

D$66 

M$66 


1 

5 
1 
1 
1 
6 
1 
6 
2 
1 
6 


100 


1 


DLOG10 


DLOG10 


F$AT 

DLOG2 

M$66 


1 
1 

1 


10 


1 


DMAX1 


DMAX1 


L$66 
H$66 
S$66 


3 
2 
1 


40 


3 


DMIN1 


DMIN1 


L$66 
H$66 
S$66 


3 

2 

1 


40 


3 



D-6 



AM74 











Approx. 




Primary- 


Entry- 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words 1Q ) 
20 


Number 


DMOD 


DMOD 


F$AT 


1 


1 






L$66 


1 










D$66 


1 










H$66 


1 










DINT 


1 










M$66 


1 










S$66 


1 










N$66 


1 






DSIGN 


DSIGN 


F$AT 

L$66 

N$66 


1 
3 

1 


20 


1 


DSIN 


DSIN 


F$AT 

DABS 

M$66 

H$66 

C$61 

C$16 

N$66 

A$66 

MOD 

L$66 

S$66 


1 
1 

9 
5 
1 
1 
3 
7 
1 
8 
2 


130 


1 


DSQRT 


DSQRT 


F$AT 

L$66 

C$62 

H$22 

SQRT 

C$26 

H$66 

D$66 

A$66 

A$81 


1 

2 
1 
1 
1 
1 
1 
1 
1 
1 


40 


1 


E$ll 


E$ll 


ARG$ 
M$ll 
F$ER 


1 
2 
1 


100 


3S 


E$11X 


E$ll 


ARG$ 


1 


110 


3H 




E$11X 


F$ER 


1 






E$21 


E$21 


ARG$ 
M$22 
D$22 


1 
1 
1 


50 


3 


E$22 


E$22 


ARG$ 
ALOG 
M$22 
EXP 


1 

1 
1 
1 


30 


3 


E$26 


E$26 


F$AT 

C$26 

H$66 

DLOG 

M$66 

DEXP 


1 
1 
2 

1 
1 

1 


30 


1 



D-7 



AM 74 











Approx. 




Primary 


Entry- 


Subroutines 


Number of 


Storage 


Tape 


Name 


Points 


Called 


References 


(Words ) 
60 


in uivioe r 


E$51 


E$51 


F$AT 


1 


2 






H$55 


3 










IABS 


1 










L$55 


4 










M$55 


1 










D$55 


1 






E$6l 


E$61 


F$AT 

H$66 

L$66 

D$66 

D$ll 

M$ll 

M$66 


1 

5 
5 
1 
2 
1 
2 


70 


1 


E$62 


E$62 


F$AT 

H$66 

DLOG 

M$62 

DEXP 


1 
2 
1 
1 
1 


30 


1 


E$66 


E$66 


F$AT 

H$66 

DLOG 

M$66 

DEXP 


1 
2 
1 
1 
1 


30 


1 


EXP 


EXP 


ARG$ 

N$22 

M$22 

S$22 

A$22 

D$22 

F$ER 


1 
2 
6 
3 
2 
2 
1 


230 


3 


F$AR 


See F$IO 








5 


F$AT 


F$AT 






58 


5 


F$B5-9 


F$B5 
F$B6 


C$MR 


6 


26 


6 




F$B7 
F$B8 
F$B9 
F$BN 










F$CB 


See F$IO 








5 


F$D5-9 


F$D2 


o$ps 


2 


34 


6 




F$D5 


0$ME 


6 








F$D6 












F$D7 












F$D8 












F$D9 












F$DN 











D-* 



AM 74 















Approx. 




Primary- 


Entry 


Subroutines 


Num 


ber 


of 


Storage 


Tape 


Name 


Points 


Called 


Refe 


rences 


(Words i() ) 

37 


Number 


F$ER 


F$ER 


AC5 




2 




6 




F$HT 














F$F5-9 


F$F5 
F$F6 
F$F7 
F$F8 
F$F9 
F$FN 


C$BR 

F$ER 




1 
1 




41 


5 


F$GA 


F$GA 


F$ER 




1 




18 


4 


F$GC 


F$GC 










14 


4 


F$HT 


See F$ER 












6 


F$IO 


F$IO 

F$CB 

F$L1 

F$L2 

F$L3 

F$L5 

F$L6 

F$AR 


F$ER 




2 




1356 


5 



F$R1 



F$R2 



F$R3 



F$R5-9 



F$Rn 



F$TR 



F$R1 



F$R2 



F$R3 



F$R5 
F$R6 
F$R7 
F$R8 
F$R9 

F$Rn 



F$TR 



F$IO 
I$AA 
I$AB 

F$IO 
I$PA 
I$PB 

F$IO 
I$CA 
I$CB 

F$IO 
I$MA 
I$MC 



F$R1 
F$R2 
F$R3 
F$R5 
F$R6 
F$R7 
F$R8 
F$R9 

F$W1 

F$AR 

F$CB 

F$L6 

AC1 

AC 2 

AC 3 



21 



21 



21 



80 



45 



198 



D-9 



Name 



F$W1 



F$W2 



F$W3 



F$W4 



F$W5-9 



F$Wn 



H$66 



IABS 
IDIM 
ID INT 
IFETCH 



Flntrv 
Points 

F$W1 



F$W2 



F$W3 



F$W4 



F$W5 
F$W6 
F$W7 
F$W8 
F$W9 

F$Wn 



FLOAT FLOAT 

H$22 H$22 

H$55 H$55 



H$66 



IABS 
IDIM 
See IFIX 
IFETCH 



Snhrontines 
Called 

F$IO 

0$AP 

0$AC 

0$AF 

0$AB 

F$IO 
0$PF 

o$pp 
o$pc 

0$PB 

F$IO 

0$CH 

0$CB 

F$IO 
0$LF 
0$LP 
0$LO 

F^jIO 
0$MC 
C$8T06 
0$MA 

F$W1 
F$W2 
F$W3 
F$W4 
F$W5 
F$W6 
F$W7 
F$W8 
F$W9 

C$12 

ARG$ 

ARG$ 
AC1 
AC2 
AC 3 
AC4 

ARG$ 
AC1 
AC2 
AC3 



Number of 
References 



Approx. 

Storage 

(Words ) 

94 



80 



39 



36 



41 



ARG$ 



20 

10 

20 

10 



Tape 
Number 



10 


4 


10 


5 


20 


3 



4 
5 



D-10 



AM 74 



Primary- 
Name 


Entry 
Points 


Subroutines 
Called 


Number of 
References 


Approx. 

Storage 
(Words iQ ) 

10 


Tape 
Number 


IFIX 


ID INT 

INT 

IFIX 


L$22 
C$21 


1 
1 




4 


INT 


See IFIX 












ISIGN 


ISIGN 








20 


5 


ISTORE 


ISTORE 


F$AT 


1 




10 


4 


L$22 


REAL 
L$22 


ARG$ 


1 




10 


5 


L$33 


L$33 








10 


5 


L$55 


L$55 


ARG$ 


1 




20 


4 



L$66 



M$55 



L$66 



LOC 


LOC 


M$ll 


M$ll 


M$11X 


M$ll 
M$11X 


M$22 


M$22 
D$22 
DIV$ 


M$22X 


M$22 


M$52 


M$52 



M$55 



AC1 
AC2 
AC 3 

AC4 

ARG$ 
AC1 
AC 2 
AC3 



ARG$ 
F$ER 

ARG$ 
F$ER 

N$22 

ARG$ 

F$ER 

F$ER 

F$AT 

H$55 

SUB$ 

L$22 

M$22 

H$22 

L$55 

F$AT 

H$55 

SUB$ 

L$22 

M$22 

H$22 

S$22 

N$22 

A$22 

L$55 



1 
1 
2 
2 
2 
2 
1 

1 
1 
10 
4 
4 
4 
1 
1 
1 
1 



20 



10 
110 

50 

330 

130 

50 



110 



4 
5S 

5H 

5S 

5H 
2 



D-ll 
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Primary- 
Name 


Entry 
Points 


Subroutines 
Called 


Number of 
References 


Approx. 

Storage 

(Words^) 

20 


Tape 
Number 


M$62 


M$62 


F$AT 
H$66 
DBLE 
M$66 


1 

1 
1 
1 


1 


M$66 


See A$66 








35 


M$66X 


See A$66X 








34 


MAXO 


A MAXO 
MAXO 


FLOAT 


1 


40 


3 


MAX1 


AMAX1 
MAX1 


L$22 
H$22 
S$22 
IFIX 


2 
2 
1 
1 


50 


3 


MINO 


AMINO 
MINO 


FLOAT 


1 


30 


3 


MINI 


AMIN1 
MINI 


L$22 
H$22 
S$22 
IFIX 


2 
2 
1 
1 


50 


3 


MOD 


MOD 


D$ll 
M$ll 


1 
1 


20 


4 


N$22 


N$22 






10 


5 


N$33 


N$33 






10 


4 


N$55 


N$55 


H$55 
SUB$ 
L$22 
N$22 
H$22 
L$55 


1 
2 
2 
2 
2 
1 


30 


2 


N$66 


N$66 


AC1 
AC2 
AC 3 


1 
1 

1 


30 


4 



OVERFL OVERFL 

REAL See L$22 

S$22 See A$22 

S$22X SeeA$22X 

S$52 S$52 



AC 5 



S$55 



S$55 



F$AT 

H$55 

L$22 

S$22 

H$22 

L$55 

F$AT 

H$55 

SUB$ 

L$22 

S$22 



4 
2 
2 



20 



30 



40 



5 

5 
5S 
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Primary 
Name 


Entry- 
Points 


Subroutines 
Called 


Number of 
References 


Approx. 
Storage 
(Words 1Q ) 


Tape 
Number 






N$22 
H$22 
L$55 


2 
2 




S$62 


S$62 


F$AT 

H$66 

DBLE 

S$66 

N$66 




20 


1 


S$66 


See A$66 








3S 


S$66X 


See A$66X 








3H 


SIGN 


SIGN 


L$22 
N$22 


2 
1 


20 


4 


SIN 


COS 

SIN 


ARG$ 

N$22 

M$22 

S$22 

A$22 


1 

2 
7 

1 
4 


190 


3 


SLITE 


SLITE 

SLITET 

SSWTCH 


ARG$ 
L$33 


3 
1 


70 


5 


SLITET 


See SLITE 










SNGL 


See C$62 










SQRT 


SQRT 


ARG$ 

DIV$ 

D$22 

A$22 

F$ER 


1 

1 

1 


70 


3S 


SQRTX 


SQRT 
SQRTX 


ARG$ 
D$22 
A$22 
F$ER 


1 

1 
1 
1 


80 


3H 


SSWTCH 


See SLITE 










SUB$ 


SUB$ 


M$ll 
F$ER 


3 
1 


130 


4 


TANH 


TANH 


L$22 
EXP 
A$22 
H$22 
D$22 


1 
1 
2 
1 
1 


60 


3 



Z$80 



Z$80 



AC1 



20 



D-13 



APPENDIX E 
ERROR MESSAGES 



Error 

Message 

AD 
BF 

DL 
DT 
DZ 
EQ 

EX 
FE 
GO 

II 



IM 

IN 
IZ 
LG 
MD 

PZ 

RI 

SA 
SD 
SM 

SQ 



Condition 

Over /underflow in double -precision 

End-of-file mark encountered while unit 
backspacing a record. 

Negative or zero argument 

Both arguments are zero 

Division by zero 

Exponential overflow adding integer to 
double -precision exponent 

Exponential overflow during exponentiation 

Format error 

Incorrect control variable in a GO TO 
statement 

First argument zero, second argument 
negative 

I > 2 and J > 15, or 

I< -2 and J > 15 

Over /underflow during integer 
mul ti pli c atio n 

Input error 

Integer division by zero or -32, 768/ -1 

Log of negative or zero argument 

Double -precision multiplication or 
division over /underflow 

Double-precision division by zero 

Integer too large when converted from 
real to integer 

Arithmetic overflow (result > 2**127) 

Divisor unnormalized 

Arithmetic overflow during multiplication 
or division 

Negative argument 



Subroutine 
A$66, S$66, A$66X, S$66X 
F$F5-9 

DLOG, DLOG10, DLOG2 

DATAN2 

D$22, D$22X 

A$81 

EXP 
F$IO 
F$GA 

E$ll, E$11X 



M$ll, M$11X 

F$AR 

D$ll, D$11X 

ALOG, ALOG10, ALOGX 

D$66, M$66, D$66X, M$66X 

D$66, D$66X 
C$21 

A$22, A$22X 

D$22 

M$22, M$22X, D$22X 

SQRT, SQRTX 
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